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Abstract 

The  transition  solution  across  a  plane  magna tohydrodynamlc 
shock  of  finite  strength  In  nonallgned  flows  In  a  perfect  gas 
is  Investigated.   A  unified  point  of  view  Is  taken  to  analyze 
shocks  In  both  unsteady  one -dimensional  flows  and  steady  two- 
dimensional  flows.   Sets  of  conversion  formulas  connecting 
pertinent  parameters  In  these  two  cases  are  given.   We  present 
a  convenient  parametric  representation  of  the  shock  solution 
from  which  an  analytical  solution  expressed  In  terms  of  any 
desirable  shock  strength  parameter  can  be  readily  obtained  by 
performing  a  straightforward  calculation.   Major  emphasis  Is 
then  placed  on  the  analytical  and  numerical  solutions  In  terms 
of  the  flow  deflection  angle  across  the  shock,  the   parameter 
most  often  prescribed  In  steady  two-dimensional  flows;  entire 
ranges  of  upstream  parameters  are  covered.   Limiting  solutions 
for  several  special  cases  are  derived.   Analyclcal  solutions 
In  terms  of  the  miagnetlc  field  deflection  angle  across  the 
shock  are  also  given.   We  Illustrate  the  connection  of  the 
present  general  solution  with  several  other  m.ajor  results  In  the 
literature.   Dynamical  stability  criteria  associated  with 
general  magnet ohydrodynamlc  shocks  are  discussed. 
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I.   Introduction 

The  transition  solution  of  the  physical  state  across  a  plane 
magnet ohydrodynamic  (MHD)  shock  In  a  compressible.  Infinitely 
electrically  conducting,  invlscid,  non-heat  conducting  perfect 
gas  in  the  presence  of  an  arbitrarily  oriented  uniform  magnetic 
field  is  investigated.   In  the  absence  of  all  dissipation 
mechanisms,  the  shock  appears  as  a  surface  of  discontinuity 
with  infinitesimal  thickness  across  which  the  physical  state 
undergoes  a  sudden  jump  according  to  conservation  laws.   Corres- 
ponding to  each  given  upstream  state  ahead  of  the  shock,  there 
exists  a  family  of  downstream  states  behind  the  shock  charac- 
terized by  a  single  shock-strength  parameter  for  which  any 
convenient  downstream  physical  quantity  (including  the  shock 
speed  or  shock  inclination)  can  be  taken.   Thus  the  transition 
solution  consists  of  determining  the  downstream  state  in  terms 
of  a  prescribed  shock-strength  parameter,  together  with  the  known 
upstream  state.   Because  of  the  presence  of  Maxwell  stresses, 
the  flow  direction  changes  across  a  MHD  shock  in  general,  even 
if  it  is  normal  to  the  shock  surface  upstream  (see  Section 
II -A)-,  hence  general  MHD  shocks  are  always  oblique. 

MHD  shocks  play  an  essential  role  in  both  unsteady  and 
steady  MHD  flow  problems  when  compressibility  effects  are  duly 
taken  into  account.   Their  solution  has  been  investigated  by 
many  authors.   The  conservation  laws  across  a  MHD  shock  were 
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were  first  formulated  by  de  Hoffmann  and  Teller.   Earlier 

2  4 
works    are  mostly  concerned  with  the  classification  of 

different  types  of  MHD  shocks,  their  general  analyticai  ifjo^o- 

perties  and  explicit  transition  solutions  for  various  simple 

special  cases.   MHD  shocks  in  aligned  flows  in  which  the 

magnetic  field  has  the  same  direction  as  that  of  velocity  on 

5-12 
both  sides  of  the  shock  were  studied  extensively  later. -^ 

The  problem  is  much  simpler  in  this  case;  detailed  basic 

properties,  their  analytical  and  numerical  solutions,  were 


1.  F.  de  Hoffmann  and  E.  Teller,  Phys .  Rev.  80,  692  (1950). 

2.  R.  Liist,  Z.  Naturforsch.  8a,  277  (1953). 

3.  R.  Liist,  Z.  Naturforsch.  10a,  125  (1955). 

4.  K.  0.  Friedrlchs,  "Nonlinear  Wave  Motion  in  Magneto -hydro- 
dynamics," Los  Alamos  Report  LAMS  2105  (195^)-   See  also 
K.  0.  Priedrichs  and  H.  Kranzer,  "Nonlinear  Wave  Motion," 
Courant  Institute  of  Mathematical  Sciences,  New  York  Univ. 
Report  NYO-6486-VIII  (1958). 

5-   H.  L.  Heifer,  Astrophys.  J.  II7,  177  (1953). 

6.  J.  Bazer  and  W.  B.  Ericson,  Astrophys.  J.  129,  758  (1959). 

7.  W.  B.  Ericson  and  J.  Bazer,  Phys.  Fluids  3,  631  (i960). 

8.  J.  Bazer  and  W.  B.  Ericson,  Proc .  of  Symp .  on  Electromag- 
netics and  Fluid  Dynamics  of  Gaseous  Plasma,  Polytechnic 
Institute  of  Brooklyn,  387  (1962). 

9.  M.  N.  Kogan,  PMVI  (Appl.  Math,  and  Mech.)  23,  92  (1959). 

10.  M.  N.  Kogan,  PMM  (Appl.  Math,  and  Mech.)  25,  784  (1959). 

11.  H.  Cabannes,  Rev.  Mod.  Phys.  32,  973  (i960). 

12.  R.  V.  Polovin,  Sov .  Phys.  Uspekhi  3,  677  (1961).   Extensive 
bibliographies  are  included  in  this  review  paper. 
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obtained.   Obviously  actual  MHD  flows  are  nonaligned  in  general 
Common  justification  for  the  study  of  only  the  restricted,  and 
therefore  simpler,  case  of  aligned  flows  states  that  even 
though  a  general  MHD  flow  is  nonaligned,  one  may  always  trans- 
form it  into  an  aligned  flow  by  translating  the  inertia  frame 
along  the  shock  surface  with  an  appropriate  uniform  speed  to 
eliminate  the  velocity  component  normal  to  the  magnetic  field 
upstream,  hence  also  that  downstream.   Strictly  speaking,  this 
is  possible  only  for  unsteady  one -dimensional  flows,  although 
needless  to  say  the  inverse  transformation  is  usually  very 
tedious.   In  steady  MHD  flows  the  relative  velocity  between 
the  oncoming  flow  and  the  source  of  disturbances  (e.g.,  the 
flow  velocity  with  respect  to  a  fixed  object  dovjnstream)  is 
prescribed  a  priori  and  it  is  invai'iant  under  any  transforma- 
tion of  inert ial  frame;  one  no  longer  possesses  the  extra 
degree  of  freedom  of  transforming  from  a  nonaligned  flovj  to  an 
aligned  flow.   Besides,  the  shock  inclination  is  not  known  and 
is  to  be  determined.   It  is  well  known  in  steady  MHD  flows 
that  if  the  flow  is  aligned  at  one  point,  it  should  be  so 
throughout  the  entire  flow.   Therefore  nonaligned  steady  flows 
(B  X  U  =  E  /^  0)  and  aligned  steady  flows  (B  x  u  =  E  -  O) 
represent  two  separate  and  distinct  physical  phenomena.   The 
solution  for  MHD  shocks  in  aligned  flows  actually  corresponds 
to  a  very  degenerate  limit  of  the  general  solution  for  MHD 
shocks  in  nonaligned  flows  (see  Section  II-B);  it  cannot 
describe  a  wide  variety  of  distinct  features  inherent  in 
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general  MHD  shocks  and  Is  by  no  means  representative.   To 
manifest  this  by  a  simple  example,  let  us  consider  shocks  of 
Infinitesimal  strength  (i.e.,  the  small  disturbance  waves  of 
characteristics)  In  steady  MHD  flows.   It  follows  immediately 
from  the  characteristic  locus  diagram  that  in  aligned  flows 
two  slow  shocks  are  always  directed  upstream  and  two  fast 
shocks  are  always  directed  downstream,  both  symmetric  with 
respect  to  the  flow  direction,  whereas  In  general  nonaligned 
flows,  downstream-directed  slow  shocks  and  upstream-directed 
fast  shocks  exist  and  the  symmetry  property  is  lost. 

Solutions  for  MHD  shocks  in  general  nonaligned  flows^~  ' 
have  so  far  received  only  relatively  little  attention,  mainly 
because  of  the  Inherent  analyti.cal  complexities  which  arise 
from  its  Involving  three  independent  parameters  characterizing 
any  given  upstream  state,  in  addition  to  the  ratio  of  specific 
heats  (i.e.,  the  adiabatic  exponent),  y,    which  characterizes 
the  gas.   As  a  result  the  explicit  solution  and  the  detailed 
fundamental  properties  of  general  I^fflD  shocks  are  only  partially 


13.   J.  D.  Cole,  in  Magnet odynamics  of  Conducting  Fluids,  ed 
D.  Bershader,  Stanford  University  Press  (1959). 
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M.  I.  Kiselev,  Sov.  Phys . ,  Doklady  ^,  517  (1959) 


15.  M.  I.  Kiselev  and  N.  I.  Kolosnltsyn,  Sov.  Phys.,  Doklady 
5,   246  (i960). 

16.  J.  E.  Anderson,  Magnet ohydrodynamlc  Shock  Waves,  MIT  Press 
(19d3). 

17.  Y.  M.  Lynn,  "Moderately  Strong  Oblique  Shocks  and  Simple 
Waves  m  Steady  Magnetohydrodynamic  Plow,"  Courant  Inst. 
?L   Mathematical  Sciences,  New  York  Univ.  Report  NYO-10436, 
MP -3b  (1964). 
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known  at  present.   Another  Important  point  to  observe  is  the 
choice  of  the  shock  strength  parameter.   In  previous  works 
this  is  mostly  dictated  only  by  analytical  simplicity  without 
due  regard  to  the  available  natural  parameter  actually  pro- 
vided by  the  boundary  condition  in  the  specific  flow  problem. 
Although  the  shock  solution  given  in  terms  of  these  analytically 
convenient  but  physically  improper  strength  parameters  do, 
admittedly,  assume  a  simpler  expression,  the  physical  picture 
is  totally  obscured.   In  order  to  describe  MHD  shocks  in  actual 
flow  problems  it  is  essential  to  obtain  the  correct  representa- 
tion of  the  shock  solution  expressed  in  term.s  of  an  appropriate 
shock  strength  parameter  in  accordance  with  the  known  boundary 
condition.   This  is  necessary  not  just  because  the  solution 
may  be  used  directly  and  the  corresponding  dynamical  stability 
criteria  can  be  readily  examined,  but  also  because  only  then 
the  structure  of  the  shock  solution  is  clearly  displayed,  and 
the  actual  physical  situation  of  the  flow  could  be  of  help  in 
selecting  a  pertinent  one  out  of  several  possible  solutions. 
Moreover,  the  range  of  the  given  physical  parameter  over  which 
a  shock  solution  exists  C8,n  be  determined  conveniently  from  this 
proper  representation.   Due  to  the  complexity  of  the  solution 
the  transf orm.ation  of  the  analytical  solution  from  one  form  to 
another  expressed  in  terms  of  a  new  shock  strength  pa.rameter 
usually  presents  a  difficult  task. 

The  equations  for  the  tra-nsition  solution  across  a  plane 
MHD  shock  in  unsteady  one -dimensional  and  steady  two-dimensional 
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flows  are  substantially  the  same  when  they  are  analyzed  in  a 
coordinate  system  fixed  to  the  shock.   Major  differences 
between  them  lie  in  the  available  upstream  parameters  and 
the  shock  position  parameter  in  accordance  with  each  of  their 
respective  physical  situations,  as  well  as  the  application  of 
the  solution  to  describe  the  physical  flow  in  each  case,  e.g., 
the  limited  region  of  influence  ahead  of  the  shock  on  each 
side  of  the  flow  velocity  vector  which  is  directly  associated 
with  the  direction  of  the  latter  enters  into  consideration  only 
in  steady  two-dimensionsl  flows  but  not  in  unsteady  one -dimen- 
sional flows.   The  dynamical  stability  criteria  for  these  two 
cases  are  also  different.   The  above  properties  have  not  been 
emphasized  previously  and  no  explicit  illustration  of  the 
connection  between  solutions  for  these  two  different  flows 
has  been  given. 

An  attempt  is  made  in  this  paper  to  remove  the  above 
mentioned  Inadequacies.   The  exact  transition  solution  for 
non-relativistic,  plane,  compressive  MHD  shocks  in  the  most 
general  case  is  studied.   In  Section  II  the  conservation  laws 
across  a  general  MHD  shock  are  formulated  and  analyzed  from 
the  point  of  view  of  that  in  steady  two-dimensional  flows. 
The  solution  is  given  in  a  convenient  parametric  representation 
from  which  the  explicit  shock  solution  expressed  in  terms  of 
any  desirable  shock  strength  parameter,  as  suggested  by  the 
specific  flow  in  consideration,  can  be  readily  obtained  by  per- 
forming a  straightforward  algebraic  elimination  process. 
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simple  conversion  formulas  connecting  pertinent  parameters  for 
shocks  In  steady  two-dimensional  flows  and  In  unsteady  one- 
dlmenslonal  flows  are  given;  they  serve  as  a  link  to  unify  the 
shock  transition  solutions  for  these  two  flows.   Major  empha- 
sis Is  then  placed  on  the  analytical  and  numerical  solutions 
In  terms  of  the  flow  deflection  angle  across  the  shock  which 
Is  most  often  prescribed  either  directly  or  Indirectly  from  the 
boundary  conditions  In  steady  two-dimensional  MHD  flows. 
Limiting  solutions  at  several  special  values  of  the  basic  par- 
ameters are  then  derived.   The  representation  of  the  shock  so- 
lution In  terms  of  the  magnetic  field  deflection  angle  across 
the  shock  Is  also  given.   Section  III  Illustrates  the  con- 
nection of  the  present  general  solution  with  several  other 
major  results  In  the  literature.   The  dynamical  stability 
criteria  associated  with  general  MHD  shocks  are  discussed  In 
Section  IV.   This  paper  Is  concluded  by  several  general 
discussions  In  Section  V. 
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II.   Formulation  and  Solution 

The  equation  of  state  for  a  perfect  gas  (sometimes  called 
polytropic  gas)  is  given  by  P  ==  A(S)p'^  where  y ,    the  ratio  of 
specific  heats  (or  the  adiabatlc  exponent),  is  constant,  P  and  p 
are  pressure  and  density  respectively,  A(S)  ^Pq/pI  exp(S-S  /C  ) 
is  a  known  function  of  entropy,  S,  with  the  subscript  "o" 
denoting  values  at  a  certain  given  reference  state, and  C   is 
the  specific  heat  of  gas  at  constant  volume.   The  state  of  the 
gas  at  any  point  in  the  flow  is  completely  described  by  local 
physical  variables  (p ,  P,  U,  B)  where  tf  and  B  are  generally 
three-dimensional  velocity  and  magnetic  field  vectors. 

The  non-relativistic  limit  of  the  well-known  de  'Hoffmann- 
Teller  relations  across  a  plane  stationary  MHD  shock  (or  the 
generalized  Rankine-Hugonlot  relations  in  magnetohydrodynamics) 
consists  of  the  following  system  of  equations  which  described 
the  conservation  of  mass,  of  normal  and  tangential  momenta,  of 
normal  magnetic  field  strength,  of  tangential  electric  field 
strength,  and  of  energy  across  the  shock  surface.  "^"^ 

[PU^]  =  0  (la) 

.2   ^2 


p       B„,  -  B 


fP^n^T  -  ~~^   -   0  (ic) 

-  8  - 


[Dj^gj  -  UjB^]  =  0  (le) 

["n[|("n  +  "?)+7^p]^H""^  ■  ("T-^T)Bn!]  =0    (If) 


Where  [   ]  denotes  the  jump  of  the  quantity  Inside  it  across 
the  shock  surface,  and  the  subscripts  n  and  T  denote  normal 
and  tangential  components  of  the  vector  with  respect  to  the 
shock  surface.   The  equation  of  state  for  a  polytropic  gas  has 
been  used  to  eliminate  the  internal  energy,  e,  i.e.,  e  =P/(y-1)p, 
in  the  energy  equation.   Since  Urn  and  B^n  are  two-dimensional 
vectors  in  general,  the  system  of  (la) -(if)  consists  of  eight 
scalar  equations.   Let  the  subscripts   I   and  F  denote  the 
physical  variables  ahead  of  the  shock  and  those  behind  the  shock, 
respectively.   (p_,  P,,  U  ,  B^)  are  known  for  any  given  upstream 
state . 

Define  e   to  be  the  unit  vector  normal  to  the  shock  surface 

n 

pointing  in  the  direction  of  mass  flux,  so  that  m  =  P-^U-^ 

=  p„U„   >  0.   For  general  magnetic  fields  with  B-^  >^  e^  ;^  0 
FF,n  in 


_i 


(the  special  case  of  B-^  x  e   =  0  is  discussed  in  IIB-l-e)  we 
\  tr-  In 

define  the  plane  containing  e  and  B.^.  to  be  the  (n,t)  plane  and 
e,  to  be  the  unit  vector  which  makes  a  positive  angle  (counter- 
clockwise) ^  with  respect  to  e   in  the  (n,t)  plane.   Also  define 

e   to  be  the  unit  vector  normal  to  the  (n,t)  plane  which, 
z 

together  with  e   and  e, ,  constitutes  a  right-hand  coordinate 
^  n      t 

system,  hence  e   =  e  x  e, .   Then  the  shock  discontinuity  is  in 
■^  '  z    n    t 

_\ 

the  (t,z)  plane.   It  is  evident  from  the  definitions  of  e^  and 


-  9  - 


e^  that  B^  •  e"^  =  B^  ^  =  0  or  B,  „  =  B^  ^  =  B^  .  e   . 

2  I     Z      I,Z  I,T      l,t      1 ,Z     Z 

Explicit  solutions  of  MHD  shocks  derived  from  the  conser- 
vation equations,  (la) -(if)  for  the  special  cases  of  contact 

discontinuities  with  U-.   =  Ut-,  ^  =  0  and  of  transverse  shocks 

I,n    P,n 

with  u|^^  =  Bp^P-Pp  and  p^  =  p^,  and  hence  U^^^  =   U^^^  . 

B-p  /m-P  ,  are  quite  simple  and  well  known.   Only  compressive 

shocks  (U,  /^   0  and  excluding  U^  ^  =  U^  ^  =  B-.  J^l^p-r.)   which 
X,n  -L^n    p,n    _L,irx 

correspond  to  fast  and  slow  shocks  are  investigated  in  this 
paper.   B^-   /  0  and  B^  ,  /  0  are  assumed  throughout  the  general 
analysis,  whereas  special  solutions  for  B,   =0  and  B.^  ,  =  0 
are  recovered  later  as  limits  of  the  general  solution.   Since 

^P,n  ^  -  ^F,n/^'^  ^"^  ^I,n  ^  -  ^I,n/^^^^  ^^  follows  from 
(ic)  and  (le)  together  with  (la)  and  (Id)  that  the  direction  of 

Brn  always  remains  the  same  across  the  shock  surface.   Therefore 

B^-e   =  B,-,  ^  =  0   or  B„  „  =  B^-,  ,   which  means  the  magnetic 

field  vectors  both  ahead  of  and  behind  the  shock  are  always  in 

the  (n,t)  plane.   However,  U^  is  a  three-dimensional  vector  and 

U  •  e   =  U-p   /  0  generally.   Since  B..  „  =  B„   =0  for  a  com- 

pressive  magnet ohydrodynamic  shock,  it  follows  from  (ic)  or  (le) 

together  with  (la)  and  (id)  that  U   is  always  conserved  across 

the  shock,  i.e.,  U    =  U„  ^.   Elsewhere  U   appears  only  in  the 

energy  equation  (if)  where,  since  pU  U   remains  constant  across 

the  shock,  the  corresponding  term  can  be  discarded.   Consequently 

any  constant  value  of  U^  can  be  superimposed  on  both  sides  of  the 

shock  without  affecting  the  transition  solution  of  other  physical 
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variables  across  a  compressive  magnet ohydrodynamlc  shock. 
Without  loss  of  generality  we  consider  In  the  following  only 
U  -  U  e   and  B.   All  vectors  are  then  In  the  (n,t)  plane  and 
the  physical  phenomenon  Is  essentially  reduced  to  a  plane  one. 
The  shock  surface  becomes  a  line  of  dlscontinllty  along  e.  In 
the  (n,t)  plane  and  It  may  be  called  the  shock  line  simply. 
For  later  convenience  we  define  u  s  iJ  -  U  e  ,  hence  u  =  U 
and  u,  =  U, . 

Apart  from  U  which  assumes  the  same  constant  value  on 
both  sides  of  the  shock  and  can  be  easily  taken  Into  account, 
the  state  of  the  gas  Is  now  completely  characterized  by  six 

scalar  quantities,  viz.,  p ,  P ,  u  ,  u ,  ,  B  ,  B,  ,  In  terms  of 

n   X/   n    u 

which  any  other  physical  variables  can  be  expressed.   Our 
problem  is  also  reduced  to  finding  the  explicit  solutions 
of  six  scalar  equations  in  the  system  (la)-(lf)  which  govern  the 
transition  of  the  physical  sta,te  across  a"  plane  MHD  shock  In 
both  unsteady  one-dimensional  flows  and  steady  tvjo-dimensional 
flows.   For  unsteady  one -dimensional  flows,  e   is  known  since 
it  can  only  be  in  the  prescribed  direction  of  spatial  variation. 
The  (n,t)  plane  is  then  determined  by  B-^  and  e  ,  e,  ,  and  e   are 
obtained  subsequently.   For  steady  two-dimensional  flows,  one 
is  given  the  ignorable  direction  along  which  the  physical 
state  does  not  change.   Let  e_,  be  a  unit  vector  in  this  ignorable 


direction;  the  planes  in  which  the  spatial  variation  may  take 
place  are  perpendicular  to  e^ .  Naturally  e  should  be  in  the 
plane  of  spatial  variation,  hence  e   I  e^,  which  means  that  e„ 

n  -^   (a  u 
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should  always  be  in  the  (t,z)  plane.   There  are  two  different 
types  of  two-dimensional  flows  in  magnetohydrodynamics.   In 
restricted  two-dimensional  flows  B^  J_  e^,  hence  B  _[_  e^  always, 
and  it  follows  immediately  that  e^^  1 1  ^"2  •   Therefore  the  (n,t) 
plane  is  perpendicular  to  the  prescribed  e^  and  it  is  known 
readily.   The  direction  of  e"^,  being  the  orientation  of  the 
shock  line  in  the  (n,t)  plane  (also  giving  the  orientation  of 
the  shock  surface  (t,z)  plane  which  is  perpendicular  to  the 
(n,t)  plane),  is  then  to  be  determined  by  a  single  shock 
Strength  parameter  and  the  given  upstream  state  of  the  gas.   In 
general  two-dimensional  flows  a  finite  magnetic  field  component 
normal  to  the  plane  of  variation  is  present  so  that  e  -B-j.  /=  0. 
Therefore  e_  x  e  ^   0 ,   which  means  that  the  angle  between  e_,  and 
e   in  the  (t,z)  plane  may  take  any  value  other  than  0  and  ir • 
The  determination  of  the  (n,t)  plane  from  the  presecibed  e^  is 
less  direct  while  additional  geometrical  param.eters  should  be 
considered. 

In  this  paper  only  the  shock  solution  in  the  (n,t)  plane 
will  be  investigated.   The  result  is  directly  applicable  to 
describe  MHD  shocks  in  restricted  two-dimensional  steady  flows 
and  unsteady  one -dimensional  flows.   The  use  of  this  solution 
to  describe  shocks  in  general  two-dimensional  flows  is  straight- 
forward, since  only  considerations  of  a  purely  geometrical 
nature  are  needed  in  addition.   However,  the  detailed  analysis 
is  rather  complicated  and  it  will  not  be  further  elaborated. 

We  choose  to  analyze  the  shock  solution  from  the  point  of 
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view  of  steady  flows  first.   The  resulting  solutions  are  then 
converted  to  the  form  suitable  for  describing  shock  solutions 
in  unsteady  flows.   To  this  end,  we  take  the  known  direction  of 
u_  as  a  fixed  reference  with  respect  to  which  all  angles  are 
measured,  counterclockwise  to  be  positive,  and  define  p  to  be 
the  shock  angle  made  hy  e  ,  0  to  be  the  angle  made  by  u^  (i.e., 
the  deflection  angle  of  the  flow  direction  across  the  shock), 
ifj   and  f^   to  be  the  angles  made  by  B,  and  B^  respectively  as 
shown  in  Fig.  1.   In  terms  of  polar  coordinates  with  the  direc- 
tion  of  u_  being  the  polar  axis,  we  have  u  =  (u  ,  Oj, 
Up  =  (Up,  0),  B^  =  (Bj,  -^j),  and  B^  -  (B^,  V'p)  where  in  each 
bracket  the  first  term  is  the  magnitude  of  the  vector  and  the 
second  term  is  the  angle  made  by  the  vector  with  u^.   It  may 
be  seen  from  Pig.  1  that  ahead  of  the  shock,  i^t  j^  =  u_sin  p, 
Ut-  ,  =:  u^cos  p,  B-r   ==  B-rSin(p-tT-)  and  B,  ,  =  Bt-cos(P--^  ), 

1,^      JL  -i-  >^  -L  J-  J-5^      -L  -L 

while  behind  the  shock,  ^p  v,  =  ■UpSin(p.-0 ) ,  ^p  |^  =  ^^003(^,-0), 

B-p,   =  B^s±n(B-if„)    and  B-^,  ,  =  B,-,cos  (p-T^„)  .   The  following  dimen- 
F,n     F        F        F,t     F        F 

sionless  variables  are  introduced  to  normalize  the  magnitudes 
of  the  downstream  physical  variables  with  respect  to  those 

upstream,  viz.,  p"  =  Pp/Pj^   P  =  ^j/Pj'   ^  ""  '^j/^I  ^^^     ^  "  ^l/^I' 
These,  together  with  the  angles  p,  0,  and  if-^,    constitute  seven 
dimensionless  parameters.   Any  one,  or  a  combination  of  them, 
should  be  prescribed  as  the  shock  strength  parameter;  the 
remaining  six  independent  variables  describe  completely  the  down- 
stream state  behind  the  shock.   The  given  upstream  state  is 
governed  by  three  independent  dimensionless  parameters  defined  by 
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Mj  =  Uj/a^,   kj  =  bj/a-j     and  i/^,   where  a  =  y  (dP/dpjg  =  KyP/P 
Is  the  gasdynamlc  sound  speed  and  b  =  B//TIp"  Is  the  Alfve'n  wave 
speed.   M  is  the  Mach  number  and  M,  measures  the  upstream  flow 
speed,   k  =  —  x (magnetic  pressure/gas  pressure),  hence  k 
measures  the  upstream  magnetic  field  strength.   t^_  characterizes 
the  upstream  magnetic  field  direction.   The  adiabatic  exponent, 
Y,  also  enters  here  as  a  parameter  which  depends  on  the  nature 
of  the  gas  considered  and  is  known  at  the  outset.   It  may  be 
remarked  here  that  the  streamlines  in  MHD  flows  cannot  be  replaced 
by  solid  boundaries  in  general,  for  the  magnetic  field  boundary 
condition  should  also  be  considered. 

Expressed  in  terms  of  the  dimensionless  quantities  defined 
above,  (la)-(lf)  become  '  respectively, 

sin  P  =  pu  sin(p  -  0)  (2a) 

^  +  M^sin^P  +-^  cos  2(p  -if^)    =   l  +  M^pu  sin^O-0)  +  ^—  cos  2{^-if^) 

(2b) 

M^sin  2p  -k^sin  2^-^^)  =.  M^  pu^sin^(p-e)  -k-j-B^sin  2^-^^) 

(2c) 
sin(p-T^^)  =  B  sln(p--^p)  (2d) 

sin  -if^  =   uB  slniif   -e)  (2e) 

"i      1 

^~  ^  Y^^^^"  ^  +  k^sin  T^^cos(p-^^)  (2f) 

2--2 
M^pu      p 

=  u[(— ^—  +  --^)sin(p-0)  +  k^B^sln(^p-0)cos(p--^p)]. 
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(a)   Solution  In  terms  of  p 

We  present  the  shock  solution  expressed  In  terms  of  the 
density  ratio,  p,  which  is  taken  as  the  prescribed  shock  strength 
parameter.   However,  a  shock  In  steady  flows  Is  most  conveniently 
characterized  by  Its  orientation,  e.g.,  the  Inclination  angle  of 
the  shock  line  with  respect  to  the  given  upstream  flow  direction, 
p,  which  determines  the  position  of  the  shock  line  and  marks 
geometrically  the  boundary  In  a  flow  field  separating  the  given 
upstream  state  and  the  unknown  state  of  gas  downstream.   Other 
downstream  physical  variables  to  be  evaluated  subsequently  can 
then  be  determined  In  their  respective  regions  In  the  flow.   For 
cases  where  p  Is  not  prescribed  It  Is  always  desirable  to  deter- 
mine p  first  from  the  prescribed  shock  strength  parameter  together 
with  the  given  upstream  parameters,  M^ ,  k  ,  and  i/'  .   To  this  end, 
a  relation  between  p  and  p  should  be  first  obtained. 

It  is  obvious,  by  the  definitions  of  e  ,  e  ,  and  p  that 

TT  >  p  >  0  always.   Instead  of  i/    ,    ii^,    9   and  p,  it  is  convenient 
—   —  11* 

to  use  the  tangent  of  these  angles  defined  by 


s  =  tan  1^  ,  j6  =--   tan  i^  ,   e  =  tan  6,      t  =  tan  p.    (3) 


However,  j6   and  e  characterize  only  the  inclination  angles  (or 
orientations)  of  the  unknown  downstream  vectors  B„  and  u„  with 
respect  to  u_;  the  senses  of  these  vectors  along  their  orienta- 
tions which  depend  on  the  directions  of  u-j.  and  B-j-  should  also  be 
found.   With  the  aid  of  (3)j  the  cosines  of  angles  are  given  by 
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COS  1f^  =   sg(V'j)(l+s^)  ^   COS  T^p  -  sg{if^){l+j6^)    S 

cos  0  =  sg(0)(l+e^)"2   ,   cos  p  =  sg(p)(l+T^)"2  ^ 

where  sg(^)  denotes  the  sign  of  cos  ^  and  equals  +1  If 
5  >  ^  >  -  f ,  while  It  equals  -1  If  tt  >  ^  >  |  or   -  |  >  ^  >  -tt. 
Expressed  in  terms  of  the  new  variables  defined  by  (3)  and  (4), 
(2a)  and  (2d)  become 


p(T-e) 
and 


^  -   sg(0)  ^J^  (5a) 


B  =  sg(^t'^)sg(^^) 


{---)n^  ,  (5b) 


1+s 


(T-J^)  /^  ■  -"^ 


respectively.   Since  u  and  p  are  scalars,  (5a)  yields  the  condi- 
tion for  sg(0)  that 

sg(0)(T-£)T  >  0  (6as) 

which  corresponds  to 

P  >  9  >  P-TT   .  (6b) 

Similarly,  since  B  is  a  scalar,  (5b)  yields  the  condition  for 
sg(^p)  that 

sg(V^^)sg(^p)(T-s)(T-j25)  >  0  (7a) 

which  corresponds  to 
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If  p  >  i/'j  >  P-TT  ,  then  p  >  "^p  >  P-TT 


If   TT-p  >^j>^       or    TT-p  >  t^  >  p  -2tt  (7t)) 

then  TT-p  >  V'p  >  P,    or  p-TT  >  -i^p  >  p-2Tr 


with   sg('i/^  )  being  known  from  the  given  upstream  state.   There - 

fore  the  senses  of  u„  and  B^  along  their  orientations  which  are 

characterized  by  sg(0)  and  sg('^„),  respectively,  are  determined 

from  (6)  and  (7).   In  fact  (6)  and  (7)  follow  directly  from 

(la)  and  (id)  respectively;  they  describe  the  simple  physical 

facts  as  one  may  observe  from  (la)  and  (id)  readily  that  u^ 

points  In  the  direction  of  the  positive  n-slde  of  the  shock  line 

as  Ut-  does  while  both  B.^  and  Bt-,  point  to  the  same  side  of  the 
I  I      F  -^ 

shock  line  which  can  be  either  the  positive  n-slde  or  the  nega- 
tive n-slde.   Thus  sg(0)  and  Bg{f^)   describing  the  senses  of 
u_  and  B„  respectively  can  be  found  very  easily.   Using  (5aj 
and  (5b)  to  ellm.lnate  u  and  B  In  (2b),  (2c),  (2e),  and  (2f) 
first  and  expressing  the  result  In  terms  of  the  new  variables 
defined  In  (3)  and  (^),  we  solve  and  get  the  following 
algebraic  equations  (see  Appendix  A): 


f(T,p)  =  k^p(p-l)M^T^(l+S^)(l+ST)^  +  []yi^T^(l+S^)  -  kjP(T-s)^]  X 

|2p(l+T^)[M^T^(l+s^)  -  kjP(T-s)^]  -  M^T^(l+s^)[Y+l-(7-l)p]  X 
[M^T^-  kjp-(l+T2)]j  =  0  (8a) 
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,2  2,.  .  _2s    ,  -  ,-,  .  ^2. 


t(p-1)[M^t  (1+s^)  +  k  ps(l+T  )(t-s)] 
^      "  M^t2(1+s'^)(t^+P)  +  k^p(l+T^)(T-s)(ps-T) 

M^t^(1+s^)[Pt(1+st)-(t-s)]  -  k  ps(l+T  )(t-s)' 


H^'9)  =  ^-J 


M^t^(1+s'^)[p(1+st)+t(t-s)]  -  k  pCI+t  )(t-s)' 


(8b) 


(8c) 


P(t,p)  -  1  + 


y(P-1)M?t2 


P(1-t 


^ 


X 


2  2 


1  + 


k  p(l+ST)^[2k  p(t-s)''-(p+1)M^t  (1+s  )] 


,2  2 


2TF 


2[M^T  (l+s"")  -  k^p(T-s)  ] 


-2/   -^     1 
u  (t,P)  =  2 

1+T 


B^(t,p)  = 


+ 


1  + 


k  t(p-1) (t-s) (1+st) 


t2 


M?T^(l+s^)-k^p(T-s)^ 


(8d) 


(8e) 


(1+t'^)(1+s^) 


2   ,  p(l+ST)[M^T2(l+s2)-k  (t-s)^] 
(t-s)   + 


M^T^(l+s2)-k^p(T-s)^ 


(8f) 


Note  that  sg(i/^-|-),  sg(-^  )  and  sg(9)  do  not  appear  in  the  above 
equations  because  only  scalar  physical  variables  and  the  orien- 
tations of  vectors  (the  tangent  of  their  angles)  are  Involved. 
The  complete  solution  for  the  downstream  state  should  Include 
sg(0)  and  sg(i/'„)  determined  from  (6)  and  (7)- 

As  far  as  (8a)-(8f)  are  concerned.  It  suffices  to  confine 

our  attention  to  the  consideration  of  the  angles  1^-^,  ii^   and  9 

o     I '   F 

within  one  period,  ir ,   which  vjill  be  taken  here  as  ranging  from 
-  S^  to  1^.   Equations  (8a)-(8f)  are  invariant  under  the  following 
transformation: 
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^I  ~^     ""^I  '   ^F  — >  -"^F  '  e  ^     -e    ,      p  ^  TT-p     (9a) 
which  corresponds  to 

s  ->  -s  ,   j6  -^    -j6   ,  e  ^  -e  ,   t  -^  -t  .    (9b) 

Thus,  without  loss  of  generality,  the  upstream  parameter  i/j   can 
be  further  restricted  to  the  range  from  0  to  ^  In  the  study  of 
general  shock  solutions.   The  solution  for  0  y_  1/     y   -  ^     can 
be  easily  obtained  from  the  former  by  the  use  of  the  transforma- 
tion specified  In  (9). 

Equation  (8a)  is  a  polynomial  of  degree  six  In  t ■ or  a" cubic 
equation  in  p.   For  any  prescribed  value  of  p,  t  can  be  deter- 
mined from  (8a)  which  yields  in  general  six  distinct  real  roots 
of  T.   Alternatively  one  may  also  solve  p  from  (8a)  which  yields 
generally  three  distinct  real  roots  provided  t  is  specified  as 
the  shock  strength  parameter.   A  selection  principle  based  on 
entropy  increase  consideration  according  to  the  second  law  of 
thermodynamics,  the  dynamical  stability  criteria,  and  the 
relevant  boundary  conditions  in  the  flow  are  needed  to  single 
out  an  appropriate  solution  to  describe  the  flow  in  question 
(see  Section  IV) . 

Unlike  the  case  in  one -dimensional  unsteady  flows,  shock 
lines  in  the  (n,t)  plane  in  steady  flows  to  not  extend  indef- 
initely in  both  directions  along  the  t-axls.   In  particular, 
at  the  foremost  end  of  the  extraneously  introduced  source  of 
disturbances,  say  0,  which  encounters  the  oncoming  flow  first 
(e.g.,  the  leading  edge  of  a  pointed  body  or  a  line  current 
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sheet)  only  semi -Infinite  attached  shock  lines  may  exist  on 
each  side  of  u  .   Let  (P)  denote  the  positive  side  of  u-j. 
belng  the  region  swept  by  vectors  making  angles  ranging  from 
0  to  TT  with  respect  to  u^  and  (N)  denote  the  negative  side  of 
u   being  the  region  swept  by  vectors  making  angles  ranging 
from  -TT  to  0  with  respect  to  u    .      On  the  (P)  side  the  shock 
line  originates  from  0  and  extends  indefinitely  along  the 
positive  t-axis,  whereas  on  the  (N)  side  the  shock  line  ori- 
ginates from  0  and  extends  indefinitely  along  the  negative 
t-axis.   Naturally  the  shock  angles  in  (P)  and  (N)  are  different 
in  general.   Moreover,  the  shocks  may  be  either  downstream- 
directed  if  1^  <  p  <  TT  in  (N)  or  0  <  p  <  J  in  (P),  or  upstream- 
directed  if  0  <  p  <  1^  in  (n)  of  J  <  P  <  TT  in  (P).   After  p 

-2       -2 
and  T  are  known,  e ,   /) ,    F ,    u    ,    and  B  can  be  easily  obtained 

from  the  auxiliary  equations  (8b)-(8f);  they  serve  to  determine 

the  downstream  state  completely. 

Note  that  (8a) -(8f)  are  all  polynomials  of  t;  more 

-  -2       -2 

explicit  expressions  of  s,  j?5,  P,  u  ,  and  B  given  in  terms  of 

p  alone  can  be  obtained  by  elimination  of  t  in  (8b)-(8f)  with 
the  aid  of  (8a).   Apart  from  their  advantages  for  analytical 
studies,  the  completely  explicit  expressions  are  not  necessary 
in  actual  numerical  computations,  since  the  value  of  t  for  any 
given  p  can  always  be  determined  first  from  (8a). 

In  unsteady  one -dimensional  flows  the  shock  is  most  con- 
veniently characterized  by  its  speed  of  propagation.   The  shock 
speed,  together  with  the  upstream  state  of  the  gas,  serves  to 
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determine  the  shock  completely.   Thus  In  a  coordinate  system 
fixed  to  the  shock  front,  one  is  given  not  only  M^,  k  and  s, 
but  also  p  (hence  t)  which  is  the  angle  made  by  the  flow  velocity 
ahead  of  the  shock  with  the  stationary  shock  line.   The  latter 
can  only  be  normal  to  the  prescribed  direction  of  spatial  varia- 
tion, i.e.,  e  .   Here  it  is  convenient  to  replace  t,  M^.,  k^  and 
s  by  the  following  four  dimensionless  parameters  corresponding  to 

ratios  of  various  upstream  velocity  components  with  respect  to  the 

-1  p 
stationary  shock  line  and  the  speed  of  sound 


^      (1+T^)2     '    ^      (1+^^)^ 


(10a) 


1 


b                k2(T-s)         b,                k2(i+Ts) 
^  =  sg(p)sg(T^  ) i-^ 2"^   ,   -s^  =  sg(p)sg(-^  ) ^-E ^-r 

^  ^     (1+T  )2(H-S^)2         ^  ^  (1+T  )2(1+S^)2 

or,    conversely, 

2  2  2  2 

u  o        u  u,  b  b,  u  b,     -  u,  b 

^       M^  ^    1      t        ,  n    ,      t        ^  n  t  t  n      /  n  hk  ^ 

^  =  ^'^I  =  —  +—'^1  =  —  +—'    ^   =  u   b     +  u^b^      (^°^^ 
taa  aa  nntt 

where  the  subscript  "I"  for  all  upstream  speeds  has  been  dropped 
for  brevity.   This  abbreviated  notation  will  be  employed  through- 
out the  rest  of  this  paper.   Note  that  sg(p)  and  sg('i^  )  do  not 
appear  in  (lOb).   Of  the  four  parameters  given  in  (10a),  only 
u,/a,  b  /a,  and  b,/a  are  given  parameters  characterizing  the  up- 


Because  of  the  absence  of  any  characteristic  length  or  char- 
acteristic time,  the  shock  solution  is  homogeneous  in  speeds 
and  only  various  speed  ratios  enter  as  basic  characterizing 
parameters . 
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stream  state.   They  are  Independent  of  any  translation  of  the 

coordinate  system  in  the  direction  of  e^.   u^  represents  the 

relative  velocity  between  the  shock  propagation  velocity  and 

the  normal  component  of  flow  velocity  upstream  In  any  coordinate 

system.   Thus  the  shock  speed  has  been  Incorporated  Into  u  . 

u  /a  should  properly  be  regarded  as  a  shock  strength  parameter; 

It  plays  the  same  role  t  does  In  steady  flows.   Incidentally,  In 

the  limit  of  shocks  of  very  weak  strength,  u  and  t  describe  the 

characteristics  In  unsteady  and  steady  flows  respectively;  they 

are  basically  different  from  all  other  downstream  parameters  In 

each  case . 

If  u  Is  not  prescribed  In  a  flow  problem.  It  Is  alvjays 

desirable  to  determine  u   first  from  the  prescribed  shock  strength 

n  '^  & 

parameter,  together  with  the  given  upstream  parameters,  u,/a, 
b^/a,  and  b,/a.   The  shock  coordinate  system  is  then  determined 
and  the  shock  propagation  speed  is  known.   Solutions  for  other 
physical  variables  follow  readily. 

Substituting  (10b)  into  (8a)-(8f)  we  get  the  following 
alternative  representation  of  the  solution  expressed  in  terms  of 
new  parameters  which  are  suitable  for  describing  shocks  in  unsteady 
one -dimensional  flows: 

f*(%rp)  =-p(p-l)u\2  +  (u^-pb^)  2pa2(u2-pb2) 

-  ^n'^7+l-(T-l)p](u^-pb^)   =  0    (11a) 

e(u^,p-)  =       y^fyf-<)^P^.\\^  (,,,) 

(u^+pu^)(u^-pb^)  +  p(p-l)u  u  b  b 


n  t  n  t 
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pu,b,  (u^-b^)  +  u^b^(u^-pb^) 


t  t  ^  n  n  ■ 


n  n '  n 


n' 


(lie) 


P  (^n'P")  =  1  +  -r^ ^  + 


Yu2(p-1)  j-     Pbt[2pb^  -  (P+l)%]  ] 


pa 
2 


^,  2    -K2^2 
2(u   -  pb  J 


u  (u^.P)  =  — ^  12  +  t^t  ^  2 ^^J 

u   p 


n' 


u   -  pb 

n   ^  n 


-2 


B  (u^,p) 


^Ib?. 


[■ 


P^(<  -  <)  2 


t^  n 


n  ■ 


b   I        u^  -  Pb^ 


(lid) 


(lie) 


(llf) 


sg(T^^),  sg(i/^„)  and  sg(0)  do  not  appear  in  the  above  equations 

whereas  their  consequences  can  be  easily  accounted  for  as  discussed 

2     - 
before.    (Ha)  is  a  cubic  equation  for  both  u   and  p.   With  p 

prescribed  as  the  shock  strength  parameter,  u  can  be  determined 

readily  from  (Ha)  and  It  yields  the  shock  propagation  speed. 

Alternatively,  If  the  shock  propagation  speed  Is  prescribed  as  the 

shock  strength  parameter  so  that  u   Is  known,  one  may  solve  p  from 

(lla) .   Other  downstream  variables  are  determined  subsequently 

from  the  auxiliary  equation  In  (llb)-(llf). 

The  second  law  of  thermodynamics  requires  that  the  entropy 

should  not  decrease  across  a  physical  shock,  I.e.,  S  =  S„/S-|-  >  1. 

Llist      has    shown   that   across   a  MHD   shock  p   >   1   If   and   only   If 

S  >  1.   Thus  we  may  examine  the  entropy  Increase  condition  for  a 

MHD  shock  of  finite  strength  by  considering  p  >  1  simply  (except 

2     2 
for  the  special  case  of  p  — >1  and  u  — >b   simultaneously,  see 
^  r   /       n     n  "^ 

Section  IIB-l-g  and  Appendix  D-g)  which  can  be  readily  accounted 
for  from  the  present  representation  of  the  shock  solution. 
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General  compressive  MHD  shocks  can  be  classified  into  fast 
and  slow  shocks  according  to 


B,  >  1  and  B,  <  1  (12) 


19 


respectively.  ^  We  have,  by  definition,  B^  =  B^  ^-/B-j-  ^  = 
tan{^-ilf    )/tan{^-i/^)   =  (t-s)  (H-tjZ5)/ (t-jZ5)  (1+ts)  where  j6   can  be 
eliminated  by  the  use  of  (8c).   Then 

M^T^(p-l)(l+s^)       (P-I)u^ 
B.  -  1  —  — p  o    p  p"  ~   p    _p   •    l-i-^J 

^  M^T^(1+S^)  -  pkj(T-s)^     U^  -  pb^ 

Since  p  >  1  for  all  MHD  shocks  of  finite  strength  satisfying  the 
entropy  Increase  condition,  (12)  and  (13)  yield 


u^  >  pb^  and  u^  <  pb^  (l^a) 

n^n        n^n  ^         ^ 


for  fast  and  slow  shocks  in  unsteady  one -dimensional  flows,  and 

M^T^(1  +  s^)  >  pk^(T-s)^  and  M^t^(1+s^)  <  pk  (t-s)^   (l^b) 

for  fast  and  slow  shocks  in  steady  two-dimensional  flows.   By 
the  use  of  (la)  and  (id),  (l^a)  becomes 

4,n   >  4,n     ^"^  4,n   <  4,n  (^^a) 


19.   This  basis  of  classification  can  be  refined  after  considering 
the  dynamical  stability  criteria.   The  refined  classification 
for  unsteady  one-dimensional  MHD  shocks  in  aligned  flows  with 
evolutionary  condition  taken  into  account  has  been  given  by 
Polovin  in  Ref.  12  (see  Section  rv  for  further  discussions). 
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and  (l^b)  becomes 

MpT2(l+/)  >  k^{'v-j6f      and  MpT^(l+/)  <  kp(T-^)^    (I5b) 

In  the  limit  of  very  weak  shocks  p— >  1,  (l4a  and  b)  or  (I5a  and  b) 
reduce  to  the  usual  definitions  for  fast  and  slow  characteristics 
(or  small  disturbance  waves)  given  by  limiting  values  of  u   and  x 

respectively. 

-   -       -2 
Equations  (lla),  (lid),  and  (llf)  show  that  p,  P,  and  B  do 

not  depend  on  u, .   This  Is  obvious  physically, for  the  scalar 

variables  (except  u  )  should  be  Invariant  under  any  Galilean 

transformation  of  Inertia  frame  along  the  shock  surface  (in  the 

t-dlrection) ,  whereas  only  e,  j6 ,   and  u  could  alter.   In  particular. 

If  the  fluid  velocity  ahead  of  the  shock  Is  perpendicular  to  the 

shock  line,  then  u,  =  0.   This  Is  the  case  often  encountered  In 

unsteady  one-dlmenslonal  flows;  e.g..  If  one  considers  the  shock 

propagating  Into  a  region  otherwise  at  rest,  the  shock  speed  Is 

then  given  by  u  =  u.   Setting  u,  =  0  and  u   =  u  in  (lib),  (lie), 

and  (lie)  yields  (see  also  Appendix  C) 

VtP^P-^)  .        P\^^'-^n)  -2         1     ^    .VtiPlill'       n^ 

e    =  — 2-2 '      ^   "  2-2       '      ■"      =   Z2"  +    I- — 2-2 ^  ^    ^ 

u    -pb  b    (u    -pb    )  p  u    -pb 

_2         2  2 

Note   that   u      =  p      (l+e    )    here.      For  MHD   shocks   of  finite    strength 

-         -  2      2 

p  >  1  or  p  —^   1  and  u  ^^  b   (see  Section  IIB-l-g);  hence  even 

if  u,  =  0,  e  /^  0  in  general  unless  B  =  0  or  B,  =  0.   Therefore, 

in  contrast  to  the  familiar  case  in  gasdynamlcs  (B  =  B,  =  O)  In 
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which  normal  shocks  may  exist  and  correspond  to  the  case  of 
u,  =  0  and  then  by  (l6)  e  =  0,  normal  MHD  shocks  (u^  =  e  =  O) 
occur  only  in  special  cases. of  either  B^  =  0  or  B^  =  0,  and 
general  MHD  shocks  are  always  oblique.   This  is  clear  physically 
because  the  presence  of  a  magnetic  field  generally  results  in 
different  tangential  momentum  and  hence  different  tangent 
velocity  (for  the  mass  flux  is  conserved  across  the  shock)  on 
two  sides  of  a  MHD  shock. 

p  is  not  the  shock  parameter  usually  prescribed  in  physical 
flow  problems.   Nevertheless  the  representation  of  the  shock 
solution  in  terms  of  p,  as  shown  in  (8a)-(8f)  or  (lla)-(llf), 
has  the  main  advantage  of  having  a  simple  analytical  form  from 
which  an  exact  solution  for  the  MHD  shock  transition  expressed 
explicitly  in  terms  of  any  other  desirable  shock  strength  param- 
eter, say  T],    can  be  obtained  readily  where  rj  can  be  one  of  t  (or 
u   for  unsteady  one -dimensional  flows),  e ,   ^ ,    V ,    u    ,    and  B  ,  or 
any  other  parameter  which  can  always  be  expressed  as  a  combina- 
tion of  them  including  p.   For  steady  flows,  apart  from  the 
trivial  case  of  rj  =  t,  the  procedure  is  to  obtain  first 
T)  =  n(T,p)  with  the  aid  of  (8a)-(8f).   Using  (8a)  to  eliminate 
p  in  n(T,p)  yields  f^(T,n)  =  0.   Then  f2(p,r|)  =  0  is  obtained 
from  (8a)  by  eliminating  t  with  the  aid  of  f  {i: ,r\)    =   0.   All 
other  physical  variables  expressed  explicitly  in  terms  of  r\   can 
be  determined  from  (8a)-(8f)  with  their  dependence  on  t  and  p 
being  eliminated  by  the  use  of  f^(T,T])  =  o  and  f^{p,r])    =   0, 
respectively.   Following  a  similar  procedure  as  shown  above. 
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one  may  also  obtain  from  (lla)-(llf)  an  exact  shock  solution 
for  unsteady  one-dimensional  flows  expressed  explicitly  In  terms 
of  any  desirable  shock  strength  parameter  other  than  p.   In  the 
analysis  for  this  case  u  /a  takes  the  place  of  t  and  the  given 
upstream  parameters  are  u,/a,  b  /a,  and  b,/a  Instead  of  M.^,  k-, 
and  if  J    (or  s) . 


(B)   Solution  In  terms  of  e 

(l)   Analytical  solution 

In  actual  steady  flow  problems,  the  boundary  conditions  are 
most  often  specified  In  such  a  way  that  the  flow  deflection 
angle  across  the  shock,  9    (hence  e  =  tan  9),    Is  prescribed  as 
the  shock  strength  parameter.   The  shock  angle,  p,  should  then 
be  first  determined  to  locate  the  shock  Inclination.   To  this 
end,  we  seek  a  relation  consisting  of  t  and  e  only.   This  Is 
obtained  by  elimination  of  p  from  (8a)  and  (8b);  then  (see 
Appendix  B  for  a  proper  procedure), 

f(T,e;M^,kj,s;Y)  =  (1+T£)^F^  +  {l+Tsf {r-e)F^  (1?) 

+  (l+Te)(T-e)S,  +  (T-e)^Fi^  =  0 

where  F-,  ,  Fp ,  F^,  and  Fu    do  not  contain  e  and  are  polynomials  of 
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T  given  by 

F^(t)    =  k^(T-s)^[(Y-l)M^T2(l+s^)   +  2(t-s)^] 

F^{x)    =  kJM^T^(l+s^){(Y-l)(l+s^)[2s(l+T^)   +    (2-y)t(1+st) 

-  ys(1+st)^J  +  M^T(l+s^)(T-s)^|kj[(Y+l)s(l+T^) 
+  (y-1)^t(1+st)]  +  2[2s(1+t2)  +  (2-y)t(1+st)]| 
+  2(Y-l)k^(l+ST)(T-s) 

F^(t)   =    (y-1)M^t^(1+s2)2[(y-1)t-s]    :•■..•:    [  .: 

+  M^T^(l+s2)jkj[2(T-s)[(Y+l)s^(l+T^)-Y(Y-l)(l+s^)'^^] 

+   y[ (t-1)^-s](1+st)^}+   2(1+s^)  [(y-1)(t-s)t^ 

+    (1+t^)[(y-1)t-s]]    +  M^t(1+s^)(t-s)    k^|2YT(l+ST)^ 

+     (y+1)(t-s)[  (Y-l)T(T-s)-Y(l+ST)]|+2k^^T[2(H-S^)(l+T^) 

+   y(1+st)^-2y(t-s)^]-(y+1)(t-s)(1+t^)|+4t(1+s^)(1+t^)I 
+   2k^(T-s)^Sk^[(Y+l)(T-s)2-Y(l+ST)2]    -   2(1+s2)(1+t2)J 

F^{x)    =    (y+1)    M^t^(1+s2)^[s-(y-1)t] 

+    MjT^(1+s^)(t-s) (k^[Y(T-s)+YT(l+S^)-2(l+ST)(s+T)] 

-  2t(1+S^)J+    M^T(l+S^)(T-s)^kJkj[2T(l+ST)-(Y+l)(T-s)] 

+  2[t-s+t(H-st)]'J    -  2k^(l+ST)(T-s)^  j 
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Equation  (I7)  is  a  polynomial  of  degree  8  In  t  or  a  cubic 
equation  In  e;  It  Is  the  most  general  non-relatlvlstlc  MHD  shock 
solution  relating  t  and  e  (or  p  and  9).      The  analytical  struc- 
ture of  the  solution  Is  Illustrated  by  Its  present  representa- 
tion while  various  limiting  solutions  derived  from  (17)  can  be 
readily  factorized  into  compact  expressions.   To  write  (1?)  in 
full  as  a  polynomial  either  in  t  or  in  e  for  the  purpose  of 
numerical  computations  is  immediate.   Elimination  of  p  in  the 
auxiliary  equations,  (8c)-(8f),  with  the  aid  of  (8b)  would  yield 
explicit  solutions  in  terms  of  x  and  e.   One  may  further  eliminate 
T  in  these  equations  by  the  use  of  (l?)  and  obtain  solutions  in 
terms  of  e  only.   Rather  tedious  calculations  are  Involved  in 
these  straightforward  algebraic  manipulations,  whereas  they  are 
not  necessary  in  actual  numerical  computations. 

We  present  in  the  following  the  shock  solutions  for  several 
important  special  cases  by  making  appropriate  limits  of  the 
parameters  in  (I7)  which  describes  MHD  shocks  in  steady  flows. 
The  nature  of  degeneracy  of  the  solution  in  each  of  these 
limits  as  viewed  from  the  general  solution  are  revealed  in  this 
process.   Distinct  features  characteristic  of  these  limits, 
especially  in  the  singular  cases  where  simultaneous  limits  are 
necessary,  are  discussed.   Some  of  these  limiting  solutions  have 
been  obtained  from  other  simpler  and  more  restricted  formulations 
In  literature.   The  agreement  of  all  these  diversified  limiting 
solutions  with  known  results  in  each  of  their  respective  cases 
serves  also  to  confirm  the  correctness  of  the  complicated 
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expression  in  (17).   Alternative  forms  of  each  of  these  special 
solutions  expressed  In  terms  of  (u^/a,  u^/a,  b^/a,  b^/a)  which 
are  suitable  to  describe  shocks  in  unsteady  one -dimensional 
flows,  as  well  as  the  derivations  of  the  basic  properties  sum- 
marized here,  are  given  in  Appendix  D. 

(a)  Aligned  field  (s  =  o).   This  special  case  of  MHD 
shocks  is  most  extensively  treated  in  the  literature.   Its 
solution  corresponds  to  setting  s  =  tan  i/'  =  0  in  (17)-   We  get 

f(T,e;M^,k^,o;Y)  =  t^[(y-1)M^  +  2]A  =  0  (l8) 


where 


A  =   A(t,£)  =  kjT(l+te)^  k^T(l+Te)  +  [ (2-y)M^+ (y-l )kj] (x-e ) 
+  (M^-kj) (t-e)^  [ (Y-l)MjT^+Ykj+2(l+T^)-(Y+l)kjT^] (1+Te) 


-   (Y+l)(M^-k^)T(T-£) 


Thus  three  finite  roots  of  t  in  the  case  of  general  magnetic 
field  now  degenerate  to  t  =  0.   Nontrlvial  solutions  of  t  are 
given  by  A  =  o  which  is  a  polynomial  of  degree  5  in  r  or  a 

cubic  equation  in  e.   Writing  A  =  0  in  full,  one  finds  that  It 

20 
IS  m  agreement  with  Cabannes'  result  in  ref.  11. 


20.   There  is  a  misprint  in  the  expression  for  c  in  (7)  of  Ref. 11, 
where  the  leading  term  for  the  coefficient  of  tan-^0  must  be 

Y-1  ^         v-i  h 

-   "2   -^  instead  of  -h^   x   (in  his  notations). 
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(b)  Very  weak  magnetic  field  (k  «  l).   Setting  k 
In  (17)  yields 

f(T,e;Mj,0,s;Y)  =  (l+s^ )^M^T^ ( x-e )^  M^[ (7-I )t-s ]+2 ( x-s )  W  =  0 

^  (19) 

where 

W  =   W(t,o)  =  [ (Y-l)M^T^+2(l+T^)](l+Te)-(7+l)M^T(T-e)  . 

Thus,  five  roots  of  x  degenerate  In  this  case.   Nontrlvlal  solu- 
tions are  given  by  W  =  0  which  agrees  with  the  well-known  oblique 

21 
shock  solution  In  supersonic  gasdynamlc  flows.     In  this  limit 

all  fluid  variables  are  governed  by  the  familiar  gasdynamlc  shock 

solution  as  If  the  very  weak  magnetic  field  Is  absent.   On  the 

other  hand,  the  magnetic  field  depends  on  the  shock  strength  and 

the  fluid  properties.   Thus  the  Influence  Is  only  "one  way."   It 

_2 
Is  also  of  Interest  to  note  that  B  >  1  always  In  this  case  V7hlch 

shows  that  the  shock  belongs  to  fast  shocks. 

(c)  Very  strong  magnetic  field  (k-  »  1).   Here  by  very 
strong  magnetic  field  we  mean  that  the  m.agnetlc  pressure  Is  much 

higher  than  the  gas  pressure,  and  hence  k_  »  1.   Eq.  (17)  is  a 

2 
quadratic  equation  in  k-;  retaining  only  terms  with  k  yields 

f(x,e;M^,k^»l,s;Y)  =  k^(x-s  )  (l+x^  )^[  (x-s)  (1+e^  ) 

+  y(t~£)(H-£s)]S  =  0  (20) 


21.   See,  e.g.,  H.  W.  Llepmann  and  A.  Roshko,  Elements  of  Gas- 
dynamlcs,  John  Wiley  and  Sons,  Inc.,  New  York  (1957):.  '9- ^7- 


31  - 


where 

S  =  S(T,e)  -  M^T(l+s^)[2T(s-e)+(Y+l)(T-s)e]  +  2 (t-s )^ (e-s) . 

Thus  six  roots  of  t  degenerate  in  this  case.   Nontrivlal 

solutions  for  the  shock  are  given  by  S  =  0  which  is  a  quadratic 

equation  for  t.   In  this  limit  B„  =  B^-,  while  the  fluid  proper- 

p     J. 

ties  depend  on  B,  only  through  the  orientation  of  the  latter. 
So  the  influence  is  also  "one  way,"  but  in  the  opposite  sense  as 
compared  with  the  case  of  very  weak  magnetic  field. 


(d)  Zero  flow  deflection  angle  (e  =  0).   Setting  e  =  o 
in  (17)  yields 

f(T,0;M-^,kj,s;Y)  -  P^  +  TP2  +  t^F^  +  T^F^  =  2Z-lZ2  =  0    (21) 


where 


Z^  =  Z^(t)  =  MjT^(l+S^)-M^T^(l+kj)(l+S^)(l+T2) 

+  k^(l+T^)(T-s)^ 

2  2,,^  2,.   ,,    ,.  ,  ,,  „^C._2,.._2 


Z2  =  Z2(e)  =   M^t''(1+s^)[s-(y-1)t]-(t-s)J2t^(1+s^) 

+  kjS(H-T^)[(Y+l)T-s]j   . 

Z-j_  =  0  is  a  biquadratic  equation  in  t;  its  four  solutions 
give  the  inclinations  of  shocks  of  zero  strength  across  which 
the  state  of  the  gas  does  not  change.   It  is  just  the  equation 
for  characteristics  which  can  also  be  derived  directly  from 
the  system  of  governing  differential  equations  for  the  flow. 
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The  characteristics  may  be  regarded  as  the  limiting  inclinations 
of  small  disturbance  waves  (or  shocks  of  infinitely  weak 
strength);  they  may  consist  of  two  fast  and  two  slow  waves,  two 
slow  waves,  or  four  slow  waves  In  general.   These  are  the  only 
solutions  for  shocks  of  zero  strength.   Solutions  for  weak  or 

moderately  strong  MHD  shocks  can  be  derived  readily  by  pertur- 

17 
batlons  about  these  solutions. 

Zp  =  0  is  another  biquadratic  equation  of  t.   Its  four 
solutions  describe,  however,  shocks  of  finite  strength  in  general, 
Since  k-  is  associated  with  t  which  is  the  highest  degree  term 
of  T  in  Zp  =  0  ^  T  — ^  00  should  be  one  of  the  solutions  as  k  — >  0  . 
This  corresponds  to  the  well-known  distinguished  limiting  solu- 
tion for  the  strong  branch  of  oblique  shocks  in  supersonic  gas- 
dynamic  flows,  i.e.,  p  -  '^'   ^t  =  0'  ^^^   e  -0-      It  describes  a 
normal  shock  and  is  independent  of  the  upstream  Mach  number, 
M^.   We  see  from  Zp  =  0  that  in  magnet ohydrodynamics  (k  /=  o  ) 
the  inclination  of  shocks  of  finite  strength  at  e  =  o  depends 
on  M-p,  k-|.,  and  s;  they  are  generally  not  normal  to  both  upstream 
and  downstream  flow  directions.   This  property  can  also  be  easily 
ascertained  physically,  since  the  presence  of  the  magnetic 
field  generally  yields  different  tangential  momentum  on  two 
sides  of  the  shock  surface. 

It  is  appropriate  at  this  point  to  classify  the  shock 

22 
adiabatics   In  the  (t,e)  plane,  or  equlvalently  in  the  (P,9) 


22.  The  term  "shock  adlabatlc"  is  used  here  to  describe  in  the  gen- 
eral sense  for  each  given  upstream  state  the  shock  solution 
relating  any  two  downstream  variables  (not  necessarily  p  and 
p,  e.g.,  T  and  e  are  taken  here)  of  which  one  may  be  regarded 
as  the  shock  strength  parameter. 
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plane  pertaining  to  any  given  M^,  k^,  s,  and  y,    into  two 
different  branches.  (Note  that  only  shock  adiabatics  with 
p  >  1  satisfying  the  entropy  condition  are  considered  here.) 
We  define  the  weak  branch  to  consist  of  the  parts  of  shock 
adiabatics  which  originate  from  each  of  the  solutions  for 
shocks  of  zero  strength  (i.e.,  the  characteristic  solutions 
as  discussed  above)  and  extend  up  to  include  their  respective 

first  local  maximum  (or  minimum)  that  ^  =  0,  or  equivalently , 

c\Q 

-r^  =  0.   Thus  for  any  given  upstream  state,  there  are  generally 

four  separate  weak  branches  describing  the  fast  and  slow 

shocks  on  (p)  and  (N)  sides  of  u,.   All  remaining  parts  of 

the  shock  adiabatics  are  defined  to  belong  to  the  strong 

branch,  for  they  describe  only  shocks  of  finite  strength. 

Only  if  the  shock  solution  lies  on  the  weak  branch  of  shock 

adiabatics  can  the  shock  strength  be  weakened  to  zero  by 

decreasing  |0l  monotonically  toward  zero  without  leaving  its 

own  branch.   Furthermore,  if  one  restricts  consideration  to 

shocks  of  the  weak  branch  only,  there  is  only  one  solution 

of  T  (or  p)  for  each  given  e  (or  Q)   describing  each  kind  of 

shock  as  mentioned  above.   The  shock  adiabatics  then  give 

unique  solutions  of  shocks  and  they  are  well  defined.   On  the 

other  hand,  if  the  shock  solution  lies  on  the  strong  branch, 

the  shock  strength  cannot  be  brought  to  zero  within  its  own 

branch.   It  should  first  pass  to  the  weak  branch  and  then 

approach  from  the  latter  to  the  solution  for  zero  strength 

shocks.   During  this  process,  double  or  multiple  values  of  t 
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(or  p)  for  a  certain  range  of  e  (or  0)  are  encountered. 
Apparently  the  weak  branches  of  shock  adlabatics  are  more 
interesting  physically;  they  probably  describe  most  of  the 
actual  MHD  shocks  in  nature  whenever  they  are  subjected  to 
normal  boundary  conditions.   Shocks  corresponding  to  the  strong 
branch  of  shock  adiabatics  may  occur  only  under  unusual  boundary 
conditions.   In  the  limit  of  very  weak  magnetic  fields,  the 
above  definitions  are  in  agreement  with  the  familiar  defini- 
tions for  weak  and  strong  branches  of  oblique  shocks  in  super- 
sonic gasdynamic  flows. 

(e)  Normal  upstream  magnetic  field  (B^  ,  =  o  >  or  s  =  -t   ) . 

B-r  X  e   =0  may  be  regarded  here  as  the  limit  of  B.^  ^  e  — >0^ 
In       js  In^^-^ 

so  that  e,  and  e   are  well  defined  as  usual.   This  case  appears 
t      z 

to  be  a  rather  special  situation  in  steady  flows  where  the 

orientation  of  the  shock  line  as  determined  by  a  certain  shock 

strength  parameter  specified  downstream  (i.e.,  e  here)  Just 

happens  to  be  perpendicular  to  the  given  upstream  magnetic  field 

direction.   It  yields  two  important  classes  of  shock  solutions; 

one  describes  pure  gasdynamic  shocks  and  the  other  describes 

complete  switch-on  shocks. 

B-p  ,  =  0  means  s  =  -t~  ;  setting  s  =  -t   in  (17)  yields 
1 ,  t 

f(T,e;Mj,kj,-T"^;Y)  =  t"5(i+t^ )^N^N2  =  0  (22) 

where 

N^  =  N^(T,e)  =  [ (y-1)M^t^+  2(1+T^)] (l+Te)-(Y+l)M^T(T-£) 
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Ng  =  N2(T,e)  =  [M^t(t-£)  -  kj(l+T^ ) (I+te) ] 

+  (t-e)^[Mjt2  -k^{l+x^)][2{l+x^)+{y-l)nl'v^  -  kj(T+l )  (1+t^  )  ] 

A  comparison  of  N.  =  0  with  W  =  0  In  (19)  shows  that  it  Is 
just  the  gasdynamlc  shock  solution.   Bj,  =   B-j-  across  this  shock. 
The  changes  of  other  physical  variables  are  determined  by  the 
familiar  solution  for  gasdynamlc  oblique  shock  as  If  B-j.  were 
absent.   Setting  t  =  -s~   In  N-,  =  0  yields 

2s(l+s^-M?) 

£  =  ^ ^ ^ P-  .  (23) 

M^[y-1+(y+1)s  ]  +2(l+s^) 

Eq.  (23)  Is  Independent  of  the  magnetic  field  strength  since  k 

does  not  enter  In  the  equation.   Thus  for  any  given  upstream 

state  with  M  ,  s  and  y   known  If  e  satisfies  N-,  =0  In  (22),  one 

of  the  shock  solutions  Is  a  pure  gasdynamlc  shock.   This  type 

of  shock  Is  sometimes  called  a  parallel  shock  In  the  sense  that 

B_  e      and  B„  e    . 

I    '  '      n  F    '  '      n 

2  -  2 

Np   =  0    corresponds   to   1+ts  — >0    and  M_  — >  pk    (l+s    )    slmul- 

2  2 

taneously.   This  solution  exists  only  If  k  >  1  and  M_  >  k^(l+s  ). 

It  represents  a  complete  switch-on  shock  behind  which  a  finite 

tangential  magnetic  field,  B„  ,  ,  appears,  and  obviously  It 

belongs  to  fast  shocks.   The  solutions  In  (8a)-(8f)  are  singular 

In  this  case,  and  an  appropriate  limiting  process  Is  necessary. 

The  magnitude  of  Bp  ^  depends  only  on  upstream  conditions  and  is 

given  by 
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2  2 


2       2     ^T  (       1  ^T         ) 

^'^  ^   k^  1+s       t     ^I        k^  1+s       J 


I  ^         ^  X  ..  J  ^ 


Setting  T  =  -s    In  Np  =  0  yields  the  following  quadratic 


equation  in  e: 


[M^s(l+se)-k^(l+s^) (s-e)]^  +  (l+se)^[M^  -  kj(l+s^)][ (y-I)M^ 

-  kj(Y+l)(l+s2)  +  2(l+s^)]  =  0  .  (25) 


Thus  for  any  given  M_,  k  ,  s  and  y,    a  complete  switch-on  MHDi 
shock  arises  if  e  satisfies  (25). 

The  conditions  specified  in  both  (23)  and  (25)  are  neces- 
sary but  not  sufficient,  since  only  one  of  several  possible 
solutions  of  shocks  corresponding  to  the  known  upstream  state 
and  e  gives  rise  to  the  specific  type  of  shock  discussed  here. 

Note  that  an  analogous  situation  occurs  if  B-p  ,  /^  0  but 
B-I-,  ^.  =  0  ;  it  corresponds  to  a  complete  switch-off  shock  and 
is  necessarily  a  slow  shock. 

(f)  Tangential  magnetic  field  (B^   =  o  or  t  =  s ) . 

J.  J  n 

Bt   =0  implies  B,-,   =  0  as  required  by  (la).   The  correspond- 

x  }T1  r  ,  n 

ing  shock  is  sometimes  called  a  perpendicular  shock  in  the 
sense  that  the  magnetic  fields  on  both  sides  of  the  shock  are 
perpendicular  to  e  .   B-j.   =0  corresponds  to  s  =  t;  then 
(17)  becomes 

f  (T,e;M^,k^,T;Y)  =  (2-y)M-^t^(1+t^  )^  (T-e)T  =  0  (26) 
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where 


T  =   T(T,e)  =  (Y-2)kjT(l+T  )(1+T£)''  (26) 

-[(Y-1)M^t^+  (7k^+2)(l+T2)](l+Te)(T-e)+(Y+l)M^T(T-e)2  . 


Nontrlvial  solutions  are  given  by  T(T,e)  =0  which  Is  a  poly- 
nomial of  degree  five  In  t  or  a  quadratic  equation  In  e. 
Setting  T  =  s  In  T  =  0\  yields  an  equation  from  which  e  can  be 
solved  In  terms  of  the  upstream  parameters.   Although  this  situa- 
tion appears  to  be  rather  special  In  steady  flows,  T  =  0  actually 
also  describes  MHD  shock  solutions  In  transverse  flows.   To  see 
this  let  us  consider  general  two-dimensional  steady  flows; 
transverse  flows  correspond  to  e^  |  |  B-j.  ,  and  e   |  |  B„  ,  where 
e_,  Is  the  unit  vector  In  the  Ignorable  (or  transverse)  direction. 
Thus  e„  =  e,  and  the  (Z,n)  plane  Is  the  plane  in  which  spatial 
variation  takes  place,  whereas  the  shock  surface  Is,  by  definition, 
represented  by  the  (t,z)  plane.   It  follows  from,  (ic)  that  at 

Bt-   =  Bt-,   =  0 ,  Ut-  rn  =  Ut-,  „;  henco  both  u,  and  U„  are  conserved 
I,n    F,n       I,T    P,T  t       Z 

across  the  shock.   Without  loss  of  generality  we  may  consider 
u.p   ;^  0,  u   ,  /  0  but  U    =  0  which  Is  just  the  case  presently 
considered,  although  in  m.ost  transverse  flows  one  is  concerned, 
instead,  with  the  situation  that  u    ^  0,  U_   ;^  0  but  u   ,  =0. 
Familiar  limiting  solutions  for  cases  of  very  weak  magnetic 
field  (k  —^0)  and  very  weak  shocks  (e  — >0)  can  be  obtained 
readily  from  T  =  0  in  (26).   For  the  case  of  very  strong  magnetic 
field  (k_  — >m),  t  =  0  reduces  to 
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(2-Y)eT^  +  2t  -  ye  =  0  .  (27) 

Eq.  (27)  is  a  quadratic  equation  of  t  whose  solution  Is 

2 
independent  of  M  .   Note  that  (27)  cannot  be  obtained  from 

(20)  by  taking  the  limit  of  t  — >  s ,  because  T  =  0  Is  a  linear 

rather  than  a  quadratic  equation  of  k  .   Hence  the  limiting 

processes  of  k  — >  0°  and  t  — >- s  taken  from  (17)  are  not 

interchangeable . 


(g)  Non-characteristic  solutions  with  p  =  1.   For  each 
given  upstream  state  with  M.^,  k  ,  and  s  being  perturbed, 
p  =  1  occurs  not  only  at  the  four  characteristic  solutions 
with  e  =  0  as  determined  by  Z-,  =0  in  (21),  but  also  at  two 
other  solutions  of  t  with  e  ^  0.      The  latter  solutions  appear 

l^T  (1+s  )  -^  k^i 


-  2  2     2  2 

in  the  simultaneous  limit  of  p  —^1  and  M^t  (l+s  )  — ^k  (t-s) 

2      2 
(i.e.,  u  — >-b  );  they  are  associated  only  with  slow  shocks 

across  which  7  ="6"=  1,  but  the  tangential  magnetic  field 

component  reverses  its  direction  and  u_  f^   u_  in  general. 

r      1 

Therefore  they  do  not  represent  shocks  of  zero  strength  and 
are  entirely  different  from  the  characteristic  solutions.   For 
this  case  t  is  determined  by 

[M^(l+s^)  -k^]T^  +  2k^sT  -  k^s^  =  0  (28a) 

which  yields  two  solutions.   We  have  also 


2t^(1+st) 


t(1-t  )  +  s(1+3t  ) 


(28b) 
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from  which  e  can  be  easily  determined  after  t  is  solved  from 
(28a).   These  two  new  solutions,  together  with  the  four 
characteristic  solutions,  provide  the  entire  six  end  points  of 
the  postions  of  shock  adiabatics  in  the  {r,e)   plane  which  des- 
cribe shocks  satisfying  the  entropy  condition  and  are  of  direct 
physical  interest.   Note  that  y   does  not  enter  in  all  these  six 
solutions,  and  hence  the  end  points  are  independent  of  the 
particular  kind  of  gas  in  consideration. 


(ii)  Numerical  results 

To  illustrate  the  implications  of  the  analytical  solution 
in  (17)^  numerical  results  are  obtained  by  the  use  of  an  IBM- 
709^  computer.   They  are  illustrated  in  Figs.  2  and  3- 

TT  >  P  >  0  always.   Also  it  suffices  to  consider  9   ranging 
from  -  ^  to  1^  since  (17)  is  periodic  in  0   with  the  period  t. 
Without  loss  of  generality  we  may  restrict  f     to  range  from  o  to 
^,  whereas  the  solution  for  0  >  "^-r  >  -  1"  can  be  readily  obtained 
from  the  former  by  the  use  of  the  transformation  specified  in  (9). 
In  Pigs.  2  and  3  the  solid  lines  describe  the  portions  of  the 
shock  adiabatic  satisfying  the  entropy  condition  S  >  1  (or  p  >  l), 
whereas  the  dashed  lines  describe  those  for  shocks  across  which 
S  <  1  (or  p  <  1). 

The  locations  of  the  fast  and  slow  characteristics  as  well  as 
the  non-characteristic  solutions  with  p  =  1,  are  indicated  by  03  , 
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03  ,  and  r  respectively.   For  any  given  •upstream  state  and 
gas  with  M-j- ,  k  ,  and  s  together  with  y  known,  the  corresponding 
shock  adlabatlc  on  the  (P,0)  plane  Is  obtained  by  solving  the 
polynomial  In  (l?)-   The  cubic  equation  In  e,  with  the  values 
of  T  being  prescribed.  Is  first  solved;  the  resulting  solu- 
tions are  subsequently  checked  against  that  obtained  by  solving 
the  polynomial  of  degree  eight  In  t,  with  the  values  of  e  being 
prescribed.   The  presentation  of  all  possible  solutions  covering 
entire  ranges  of  all  four  Independent  parameters  Is  obviously 
out  of  the  question.   An  attempt  Is  made  here  to  Illustrate  the 
essential  features  of  the  shock  adlabatlcs  at  certain  repre- 
sentative values  of  these  parameters.   One  may  then  Infer  from 
these  solutions  the  behavior  of  the  shock  adlabatlcs  at  other 
values  of  the  governing  parameters. 

Y  may  take  values  only  within  a  very  limited  range;  Its 
variation  does  not  alter  significantly  the  qualitative  nature 
of  the  solution  In  general.  Moreover,  all  end  points  of  the 
shock  adlabatlcs  satisfying  the  entropy  condition  (i.e.,  the 
solid  lines  In  Pigs.  2  and  3)  are.  In  fact.  Independent  of  y 
as  seen  from  Z-,  =  0  In  (21)  and  (28).  We  consider  y  -  ^  des- 
cribing monatomlc  gases  In  all  the  numerical  computations. 

Figs.  2a-f  describe  the  shock  adlabatlcs  at  M_  =  3-0  and 
1^  =  30°  with  only  k  being  different  In  various  graphs.   Thus 
the  variation  of  the  solution  due  to  the  change  of  the  strength 
of  the  magnetic  field  with  Its  direction  being  fixed  at  a 
finite  angle  with  respect  to  the  flow  direction  Is  Illustrated. 
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Guided  by  the  familiar  gasdynamic  solution  at  k^  =  o  for  which 
the  shock  adiabatic  is  antisymmetric  with  respect  to  p  =  2" 
and  is  independent  of  1/    ,   we  may  observe  from  Fig.  2a  the 
initial  emergence  of  slow  shocks  in  the  presence  of  a  weak 
magnetic  field  at  k  =  0.5-   Pigs.  2b-f  show  the  shock  adia- 
batics  as  k  increases  gradually  from  3. 0  to  ^0.0.   We  obtain 
the  equation  for  characteristic  locus  -"^  that  M-j.  =  3-0  and 
■^     =   30°;  the  upstream  flow  field  is  either  hyperbolic  with 
the  existence  of  two  fast  and  two  slow  characteristics,  or 
"hyperliptic"  with  the  existence  of  only  two  slow  character- 
istics, depending  on  whether  k   is  greater  or  less  than  8.77 
respectively.   It  is  clearly  manifested  in  Pigs.  2d  and  2e 
that  the  fast  characteristics  dieappear  as  k  increases  from 
8.0  to  10.0.   Pigs.  2e  and  2f  show  that  the  qualitative  behavior 
of  the  shock  adiabatic  remains  unaltered  as  k   increases  from 
10.0  even  up  to  ^0.0  for  the  case  of  very  strong  magnetic  field. 

For  any  fixed  k  and  any  fixed  sufficiently  small  i/    ,    let 
us  first  consider  the  case  that  M_  is  so  large  that  the  upstream 
state  is  hyperbolic;  that  is,  there  exist  two  fast  and  two  slow 
characteristics  (i.e.,  if  -u   is  outside  of  the  characteristic 
locus).   As  M_  decreases  continuously  from  that  value,  it  can 
be  seen  from  the  characteristic  locus  pertaining  to  k  that  the 
upstream  flow  field  first  becomes  "hyperliptic"  with  the  exis- 
tence of  two  slow  characteristics  only;  it  subsequently  becomes 
hyperbolic  again  with  the  existence  of  four  slow  characteristics. 


23.   Y.  M.  Lynn,  Phys.  Fluids  5,  626  (I962) 
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As  M  decreases  further  toward  zero,  the  flow  field  becomes 
"hyperllptlc"  again  with  the  existence  of  two  slow  character- 
istics only.   Naturally  it  would  be  very  interesting  to 
Investigate  here  the  variation  of  the  shock  adiabatics  as 
the  upstream  state  traverses  through  each  of  the  above 
mentioned  flow  regimes.   Since  the  inner  cusped  characteristic 
locus  corresponding  to  slow  waves  is  usually  situated  near 
ii     -  0 ,    the  upstream  state  of  flow  at  which  there  exist  four 
slow  characteristics  occurs  only  when  i/     is  sufficiently  small. 
It  is  also  of  interest  to  Investigate  the  degenerate  case  of 
i/     =   0   describing  aligned  flows  for  which  the  shock  adiabatics 
should  be  antisymmetric  with  respect  to  p  =  p-  .   Moreover, 
one  would  like  to  observe  the  gradual  approach  of  the  shock 
adiabatics  to  the  degenerate  aligned  field  solutions  as  i/j 
decreases  continuously  from  finite  values  to  zero.   The  shock 
adiabatics  for  all  cases  enumerated  above  are  Illustrated  in 
Figs.  3-   Fig.  3a  describes  the  positions  of  the  terminals  of 
-u_  relative  to  its  corresponding  characteristic  locus  for 
each  of  the  cases  considered  in  Figs.  3b-3n.   Since  the  general 
character  of  the  solutions  is  fundamentally  different  at  k  >  1 
and  k  <  1  when  if      is  small,  the  shock  adiabatics  at  k  =  0.5 
and  2.0  are  considered,  and  they  are  shown  in  Figs.  3t)-3g  and 
Pigs.  3h-3n  respectively. 

It  is  shown  in  Figs.  3f  and  3^  and  also  in  Figs.  3g  and  3n 
that  in  aligned  field  flows  there  exists  only  one  type  of 
shocks  always;  it  may  consist  of  either  two  downstream-facing 
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fast  shocks  or  two  upstream-facing  slow  shocks,  always  symmet- 
ric with  respect  to  the  flow  direction.   An  anomalous  situation 
occurs  when  M^,  k  and  ^^  are  (1.20,  2.0,  0°)  respectively. 
The  corresponding  terminal  of  -u-j.  is  at  m  as  shown  in  Fig.  3a, 
and  hence  the  flow  field  is  "hyperllptlc"  and  no  real  character- 
istics exist  except  two  degenerate  ones  coinciding  with  the 
streamlines.   Nevertheless  the  shock  solutions  satisfying  the 
entropy-increase  condition  do  exist  as  shown  in  Fig.  3m  on  which 
p  is  actually  found  to  be  always  bounded  away  from  1.   Clearly 
these  anomalous  shock  solutions  cannot  be  evolved  continuously 
from  that  for  zero-strength  shocks  (i.e.,  characteristics)  with 
9=0  and  p  =  1.   Perhaps  they  describe  physically  unrealistic 
non-evolutionary  shocks  to  be  excluded  after  the  evolutionary 
conditions  are  taken  into  account  (see  further  discussions  in 
Section  IV).   On  the  other  hand,  no  real  shock  adlabatics  satis- 
fying the  entropy-increase  condition  exist  when  the  upstream 
state  belongs  to  other  aligned  flow  reglm.es,  e.g.,  when  M.^-,  k  , 
and  1/     are  (O.85,  O.5,  0°),  (O.^O,  0.5.  0°),  and  (0.4o,  2.0,  0°) 
respectively.   Figs.  2  and  3  show  that  all  shock  adlabatics 
with  p  >  1  satisfy  (6b)  automatically.   The  portion  which  does 
not  satisfy  (6b)  actually  corresponds  to  p^  <  0  (or  p  <  0)  . 
which  is  obviously  unphysical  and  has  been  excluded  from  our 
previous  derivation  of  (6b)  from  (5a). 
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(C)   Solution  In  terms  of  e^  (=  ta.n{i/^-i'    )) 

MHD  shocks  may  be  generated  not  only  mechanically  but  also 
by  electromagnetic  means,  for  example,  by  a  magnetic  dlpole  or 
a  current  sheet  situated  downstream.   Thus  the  deflection  angle 
of  the  magnetic  field  vector  across  the  shock,  0„  =  ii^-ii^,    is 

ri      r    1 

sometimes  prescribed  by  the  boundary  condition  in  flow  problems 

as  a  natural  shock  strength  parameter. 

To  find  the  solution  in  terms  of  0_,  we  first  find  the 

ri 

solution  relating  t  and  fi   which  can  be  obtained  by  eliminating 
p  from  (8a)  and  (8c).   Thus  we  get 

g(T,;25;Mj,k^,s;Y)  =  {l-^i^)^Q^+  {l+^^f  {^ -j6)G^ 

+    (1+t/5)(t-/5)2g^+  (t-^)\  =  0  ,  (29) 

where  G-,  ,  Gp ,  G^,  and  Gh  do  not  contain  j6   and  are  polynomials 
of  T  given  by 

G-j_(t)  =  k^(T-s)5 

G^{t)    =   k^(T-s)^(l+ST)[ (2-T)M^T^(l+s^)+(Y-l)k^(T-s)^] 

G^(t)  =  (T-s)[M^T^(l+s2)-kj(T-s)2][(l+s2)[(Y-l)M^T^+2(l+T2)] 

2  2  { 

+  yk  (1+st)  -(Y+l)k  (t-s)  | 

Gi^(T)  -  -(Y+l)(l+ST)[M^T^(l+s^)-kj(T-s)2]^  . 

e„+s 
Define  e   =  tan  9^  =  tan(-^p-i^);  then  j6  =   tan  i/^  -—; 


V 
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Substitution  of  this  Into  (29)  yields 

gg(T,eg;Mj,k^,s;Y)  (30) 

3  2 

=  [1+ST+e  (t-s)  ]^G-j_  +  [l+ST+eg(T-s)]  [  T-s-eg(l+ST)  jGg 

2  3 

+  [H-ST+e  (T-s)][T-s-e^(l+ST)]  G^+[  T-s-eg(l+ST)  J^G^^  =  0 

This  Is  a  polynomial  of  degree  eight  In  t  and  a  cubic  equation 

for  e_  from  which  the  root  of  t  can  be  determined  for  each  pre- 

scribed  e  .   For  flows  aligned  upstream  s  =  0  and  (30)  becomes 

the  same  equation  as  A  =  0  In  (18),  provided  e„  Is  replaced  by 

e;  hence  0=0   and  It  Indicates  the  well-known  fact  that  the 
u 

flow  remains  aligned  downstream. 


Ill-   Connections  with  Other  Existing  Solutions 

Various  explicit  analytical  solutions  for  MHD  shocks  have 
appeared  In  the  literature.   In  this  section  the  connections  of 
some  of  the  other  major  results  with  our  general  solution  are 
Illustrated  and  discussed.   To  facilitate  direct  comparisons 
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the  notations  used  in  other  works  are  superimposed  with  a 
tilde  whenever  they  are  quoted  below. 

(a)   In  arbitrary  magnetic  field 

(l)   In  Cole's  paper   a  one -dimensional  MHD  shock  pro- 
pagating with  a  uniform  velocity  c  into  an  ideal  gas  other- 
wise at  rest  is  considered  to  be  generated  by  the  Impulsive 
motion  of  a  piston  with  a  uniform  normal  velocity  u„  into 
the  gas  behind.   In  a  coordinate  system  fixed  to  the  shock 
front,  this  corresponds  to  the  special  case  of  normal  inci- 
dence, u,  =  0 ,  discussed  in  Section  II  and  Appendix  C.   The 
correspondence  of  notations  is  u   =  c,  a  =  a^,  b   =  "^-v-' 

^t  =  ^j'    "n,F  =  ^S'  ^,P  ^  S^  ^°  ^^^^ 

v_       -■Ui-.b  b 

p  _  JL  -      F  X  y 

^  ~       ~  ~  ~    ~2  ~~   ~2 

c-Up  (c-Up)(c  -cup-b^) 

Then   1=0    in  Eq.    (Cl)    becomes 


C3_C2  =  0  (51) 


where 


0^   =    (c2-b2)(c-Up)2[2Up(B^-a2)-(Y+l)3(c2-b2)] 

-  Upby^S[Y3^-(Y+2)upH2   +    (2u|-Yb^)3  +    (Y+l)Upb^] 

2    ~2  2  2    2 

Cg  =    [2(c    -aj)-(Y^l)UpC](c    -u^c-b^) 

-  cb|[2S^-(Y+2)Up32   +    (7^2,2^2^-    ^    (^+i)up!)2]       . 


-   ^7    - 


Here  u^   Is  considered  as  the  prescribed  shock  strength 
r 

parameter,  and  the  shock  speed  c  can  be  determined  from  (31) 
C-,  =  0  gives  G  (c-b  )   =  ®  at  u„  =  0;.  it  does  not  describe 

the  compressive  shock  solution  and  must  be  discarded,   c^  = 

24 
agrees  precisely  with  Cole's  solution. 


(11)   Klselev  has  obtained  the  solution  for  general  MHD 

shocks  as  a  cubic  equation  for  B,  =  Bt-,  ^./B-^  ,  ,    together  with 

U  r  ,V        1  J  t 

u  and  the  upstream  parameters.   We  solve  for  p  from  (13)  and 


get 


u^B, 
n  t 


^^  +  b^(B^--l) 


(32) 


n   n^  t 

Substitution  of  (32)  into  (lla)  yields 

2  2         2        2    2 

-^  u^   o      u    o   2a  +(y-1)u  +Yb,   , -, 

Bj+[l-(2-Y)(l--^)]B2+(l-^)(b2 "     )^^\ 

^n  '^n  'V+1        ^t 

2        2 
K  u   2 

-  (T+1)  ^  (1  -  ^)   =  0  (33) 

\  \ 

which  is  in  agreement  with  Klselev 's  solution^  with  the 


24.  Cf.  Eq.  (2)-(7)  in  ref.  13. 

25.  Cf.  (6)  in  ref.  l4  or  (l)  in  ref.  15.   Unfortunately,  both 
of  these  equations  contain  a  misprint  in  their  last  term 
for  which  the  correct  expression  in  terms  of  original 
notations  should  read 

-(k+l)V^^h^yP^(l-^)^  . 

^Ix 
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correspondence  of  notations  being  y  ^  k,   xi     =  -u,  ,  b   =  V-,  , 

^  '     '   n    Ix   n    Ix 

The  sol-ution  for  aligned  flows  may  then  be  easily  obtained 
from  (33). 


(ill)   An  equivalent  solution  to  (lla)  was  obtained  from 

6 
a  different  approach  by  Anderson.   Dividing  (lla)  by  b   and 

expressing  It  In  terms  of  the  notations,  u  /b  =  M.-,,  b,/b  ==  M  .,  , 
^  '   n  n    Al   r  n    yl 

~        _    ~  of, 

a/b   =  M^-i  ,  and  p  =  p,  one  gets  precisely  Anderson's  solution, 
n    a  J- 

However,  no  explicit  solutions  for  other  downstream  variables 
are  available  there. 


(B)   Aligned  flows 

The  special  solution  for  aligned  flows  has  been  Inves- 
tigated by  many  authors.   The  equation  for  p  In  (lla)  now  takes 
a  simple  form  as  shown  in  (D2 ) .   Analytical  solution  for  this 
case  was  first  derived  by  Heifer.-^  A  different  and  very  com- 
prehensive study  was  later  given  by  Bazer  and  Erlcson.    The 
explicit  solution  expressed  in  terms  of  p  and  9    suitable  to 
describe  shocks  in  steady  flows  was  presented  by  Cabannes. 
Its  connection  with  our  general  solution  has  been  discussed  in 
IIB-l-a. 

(i)   The  correspondence  of  Heifer's  parameters  with  ours 


26.   Cf.  (2.^3)  on  p.  21  in  ref.  l6. 
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is   Q  =    (YK/2)(u^/a^),      s   =  u^/%.      A  -   1    -   1/kp      and 

X    =   Sp/s     =    ((A-1)p+1)/A;    one   may   solve   from  above  and  get 

p    -    (1-Ax)(l-A)    ,    K  =   1/1-Ax    ,    u^/a^   =    (2Q(1-A^  )  )/y(1+s^  ) 

and  u^/a^  =  (2Qs^ (l-Ax) )/y(1+s^ ) .   Substitution  of  these 

27 
Into  (D2)  yields  Heifer's  solution   which  Is  a  quadratic 

equation  of  A  expressed  In  terms  of  the  parameters  Q,  s, 

and  X  .   Here  the  reduclblllty  of  the  cubic  equation  of  p 

to  a  quadratic  equation  of  A  is  essentially  due  to  the  fact 

that  X  is  regarded  as  a  prescribed  parameter,  so  that 

^Z\    =  1-  X  is  then  given  and  p-x  is  Inversely  proportional 
Kp  -1 

to  1-A. 


(11)   The  correspondence  of  the  parameters  used  by 

Bazar  and  Erlcson  with  ours  is  cos  co  =  cos  p  =  "t^"  ' 

~2222~  ~  ~  ~_ 

s      =   a  /b      =  a  /ku    ,      h  =   e/(sln  co-   £   cos  cu)    and  f^   =  p-1 

The  equation  for  e  in  (D3)  becomes 


T]    sin  CDCos  o) 

1  +  T)     COS    00    -K('n+l) 


(34) 


Eliminating  e  in  the  equation  for  h  by  the  use  of  (3^)?  one 
solves  for  k  as 

K  ^  \-^^   ^°^  ^  .  •   (55) 

h(n  +  i) 

Substltution  of  {^5),    a  /u  =  ks   =  (s  (h-ii  cos  co  )  )/h(ri+l ) , 


27.   Cf.  (59)  in  ref.  5- 
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u  /u  =  sin  CD,  and  u,/u  =  cos  co  into  (D2)  yields  a  quadratic 

1-1  L/ 

equation  of  f]  from  which  t\   can  be  solved  explicitly  In  terms 

~   ~       ».  2( 

of  h,  s  ,  and  o).   This  Is  the  main  result  of  Bazer  and  Erlcson. 

Other  solutions  expressed  In  terms  of  h,  s  ,  and  co  follow 

directly.   The  crucial  step  of  the  reduction  of  a  cubic  to  a 

quadratic  equation  of  p  (=  ^-l)    lies  in  the  fact  tjiat 

K£zl  (=  £(n+llzi)  =    _  ££^   and  it  is  known  from  the  pre- 
p-1       n  h 

scribed  parameters.   Then  (D2)  becomes  a  quadratic  equation 

of  p  (or  T])  since  Kp  depends  linearly  on  p  (or  t\)   while  all 

other  parameters  are  given. 


rv.   Dynamical  Stability  Criteria 

In  previous  sections  we  were  concerned  with  the  transi- 
tion solution  across  a  general  MHD  shock  with  the  entropy 
condition  duly  taken  into  account.   The  fulfillment  of  the 
entropy  condition  alone  is  necessary  but  not  sufficient  to 
warrant  a  physically  realistic  shock.   Any  realistic  shock 
should  also  be  stable  against  small  disturbances.   These  dynam- 
ical stability  criteria  were  overlooked  in  a  number  of  works 


28.   Cf.  (19)  In  ref.  6  or  (9a)  in  ref.  8. 
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on  MHD  shocks,  perhaps  largely  because  their  counterpart  In 
gasdynamlcs  falls  to  Introduce  any  additional  restriction  on 
the  shock  transition  solution  other  than  that  already  required 
by  the  entropy  condition.   For  MHD  shocks  It  has  been  shown 
that  in  the  special  case  of  aligned  flows  the  consideration 
of  dynamical  stability  criteria  contributes  to  exclude  a  sub- 
stantial portion  of  shock  adiabatics  satisfying  the  entropy 

12 
condition.    The  corresponding  criteria  for  MHD  shocks  in 

general  nonaligned  flows  have  not  been  treated  up  to  the 
present.   In  view  of  their  underlying  importance  we  present 
in  the  following  a  formulation  of  the  dynamical  stability 
criteria  for  general  MHD  shocks  and  illustrate  the  present 
status  of  knowledge  of  these  criteria,  as  well  as  their  essen- 
tial inadequacies.   However,  no  analytical  results  are  attempted 
here,  for  the  general  problem  itself  is  complex  enough  to  merit 
a  separate  and  extensive  investigation. 

Dynamical  stability  criteria  are  usually  considered  to 
consist  of  two  types.   Firstly,  any  solution  for  a  realistic 
shock  should  possess  a  neighborhood  solution  such  that  a  small 
perturbation  of  the  shock  arising  from  arbitrary  disturbances 
on  it  should  result  in  only  a  slight  change  of  the  shock  solu- 
tion.  Mathematically  this  corresponds  to  the  requirement  of  the 
continuity  of  the  solution  of  a  realistic  shock  in  response  to 
a  continuous  change  of  boundary  data;  the  existence  of  a  slightly 
perturbed  neighborhood  solution  is  then  guaranteed.   This  condi- 
tion is  often  termed  an  evolutionary  condition  or  stability 
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against  break  up.   Secondly,  In  addition  to  the  fulfillment 
of  the  evolutionary  condition,  a  realistic  shock  should  also 
be  stable  against  time -dependent  small  disturbances  so  that 
the  temporal  development  of  the  amplitudes  of  the  perturbed 
physical  quantities  should  not  grow  indefinitely  as  time 
increases.   Both  of  these  criteria  can  be  analyzed  in  general 
terms  for  a  shock  in  an  unbounded  medium  without  referring 
to  other  local  boundary  conditions.   It  is  usually  an  easy 
matter  to  determ.ine  the  evolutionary  condition,  since  it 
demands  only  the  correct  formulation  of  a  well-posed  algebraic 
system  of  equations  governing  the  small  perturbations  of  physi- 
cal variables  pertaining  to  MHD  shocks.   To  this  end,  a  simple 
counting  of  the  num.ber  of  possible  diverging  waves  propagating 
away  from  the  shock  surface  for  all  possible  flow  situations, 
followed  by  an  examination  of  whether  this  number  being  equal 
to  the  number  of  necessary  boundary  conditions  would  suffice. 
It  should  be  emphasized  that  the  evolutionary  condition  for 
MHD  shocks  in  unsteady  one-dimensional  flovifs  is  not  to  be 
confused  with  that  in  steady  two-dimensions.1  flows.   The 
characteristics  (i.e.,  the  small  disturbance  waves)  in  steady 
two-dimensional  flows  do  not  propagate  always  in  a  preferred 
fixed  direction  in  space  as  they  do  in  one -dimensional  unsteady 
flows;  consequently  the  envelope  of  the  fronts  of  various 
normal  wave  speeds,  rather  than  the  fronts  themselves,  plays 
a  central  role  of  classifying  different  types  of  flow  situations 
In  this  respect,  the  characteristic  locus,  '    rather  than  the 
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normal  wave  speeds,  should  be  considered  in  Investigating 
the  evolutionary  condition  in  steady  two-dimensional  flows. 
On  the  other  hand,  the  determination  of  stability  against 
time -dependent  disturbances  requires  the  examination  of  the 
sign  of  the  temporal  growth  rate  of  various  amplitudes  for 
which  the  system  of  governing  differential  equations  describing 
the  flow  has  to  be  analyzed. 

Schematically,  the  entire  shock  solution  can  be  summarized 
as  follows.   Let  us  denote,  for  MHD  shocks  in  either  one- 
dimensional  unsteady  flow  or  two-dimensional  steady  flow  with 
prescribed  y   and  upstream  parameters: 
S.  =  all  transition  solutions  of  MHD  shocks  satisfying  the 

conservation  laws  described  by  (la) -(if). 
S  =  the  closed  subset  of  S.  which  satisfies  the  entropy  con- 
dition, S  >  1  (or  p  >  1). 
Sg  =  the  closed  subset  of  S.  satisfying  the  evolutionary 

condition. 
S^   =   the  closed  subset  of  S.  which  is  stable  against  arbitrary 

time -dependent  small  disturbances. 
Sj^  =  all  solutions  which  describe  physically  realistic  MHD 

shocks . 
Clearly  S„  is  our  main  concern  and  is  given  by 

2r  =  SgOSgHs^  ;  (36) 

only  S-^  and  S^  are  studied  in  this  paper.   Analytical  solutions 
of  S^  are  presented  in  Section  II.   In  Figs.  2  and  3  the  solid 
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lines  describe  S  ,  whereas  the  complement  of  S   relative 

s  s 

to  S„,  S  =  S'-S  ,  is  described  by  the  dashed  lines. 
A   s    A   s 

Evolutionary  conditions  for  MHD  shocks  in  one -dimensional 
unsteady  flows  with  arbitrarily  oriented  magnetic  fields  were 
first  considered  in  refs.  29,    50,  and  31  for  normal  and  oblique 
incident  wave  trains,  respectively.   The  results  obtained  there 

were  subsequently  applied  to  the  analytical  solution  of  one- 

12  32 
dimensional  MHD  shocks  in  aligned  flows.   '     In  this  way  the 

conclusion  that  the  evolutionary  condition  imposes  one  extra 

restriction,  viz.,  B,  >  0  (applicable  to  slow  shocks  only),  in 

addition  to  the  entropy  condition  is  valid  only  for  shocks  in 

the  degenerate  case  of  aligned  flows.   In  fact,  many  elegant 

analytical  properties  derived  for  MHD  shocks  in  aligned 

6  7 
flows    no  longer  hold  in  general  nonaligned  flows.   Naturally 

the  consequences  of  the  evolutionary  condition  for  the  solutions 

of  shocks  in  these  two  flows  are  different  in  general.   We  may 

note  that  the  evolutionary  condition  for  MHD  shocks  in  steady 

flows  has  not  been  treated  so  far. 

It  is  obvious  that  S^  3^  S   always,  because  S^  presupposes 


29.  A.  I.  Akhlezer,  G.  J.  Liubarskli,  and  R.  V.  Polovln,  Sov. 
Phys.,  JETP  8,  507  (1959). 

30.  V.  M.  Kontorovlch,  Sov.  Phys.  JETP  8  85I  (1959). 

31.  S.  I.  Syrovatskli,  Sov.  Phys.,  JETP  8,  1024  (1959)- 

32.  A.  Jeffrey  and  T.  Tanlutl,  Non-linear  Wave  Propagation, 
(Academic  Press,  196k),    p.  2T5 
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the  existence  of  a  slightly  perturbed  shock  solution  which 

belongs  to  S„.   Then  (36)  becomes 
ill 

Sj^  =  SgOs^  .  (37) 


Hence,  so  long  as  S„  alone  Is  desirable.  It  suffices  to  find 
S„  In  addition  to  S  .   The  special  case  of  S„  for  a  normal 

IS  X 

Incident  Alfven  wave  on  a  one-dlmenslonal  MHD  shock  was 

16 
given  by  Anderson.     S^  for  general  MHD  shocks  with  B^^  =  0 

33 

and  B,  =  0,    respectively,  was  studied  by  Gardner  and  Kruskal. 

The  general  case  of  B  /  0  and  B,  7^  0  Is  still  an  open  problem. 

Clearly,  the  dynamical  stability  criteria  serve  further 
to  exclude  certain  parts  of  the  shock  adlabatlcs  described 
by  solid  lines  In  Figs.  2  and  3;  presumbaly  most  of  the 
mysteriously  behaving  portions  could  be  discarded.   New 
analytical  properties  of  MHD  shocks  In  nonallgned  flows  must 
also  be  Investigated.   They  could  then  be  used  together  with 
the  deduced  results  according  to  dynamical  stability  criteria 
to  provide  several  conclusive  and  fundamental  properties  of 
general  MHD  shocks.   These  would  enable  a  substantial  refine- 
ment of  the  somewhat  crude  basis  for  the  classification  of 
fast  and  slow  shocks  considered  In  Section  II.   The  conditions 
for  the  connection  of  various  types  of  MHD  shocks  could  then 
be  derived. 


33-   C.  S.  Gardner  and  M.  D.  Kruskal,  Phys .  Fluids  7,  7OO  (196^) 
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V.   Discussions 

Only  T  (or  u  )  as  a  function  of  the  shock  strength 

parameter  Is  presented  In  detail  In  Section  II  and  Appendix 

D;  all  other  downstream  variables  follow  readily  by  the  use 

of  auxiliary  equations.   In  steady  flows  the  determination 

of  the  maximum  deflection  angle  of  streamline,  9        ,    across 

^  max 

each  kind  of  MHD  shocks  as  a  function  of  the  given  upstream 
parameters  and  y   yields  directly  the  criteria  for  the  detach- 
ment of  MHD  shocks.   In  principle,  e^„^  =  tan  9^^^  can  be 

max        max 

obtained  simply  by  eliminating  t  from  f(T,e)  =  o  and  f  =  0 

(since  de/dx  =  -  f  /f  =  o).   Several  solutions  of  9 
^  '  V     e  max 

exist  In  general,  only  a  few  of  which  pertain  to  physically 
realistic  MHD  shocks  (fast  or  slow).   The  analytical  deriva- 
tions of  9  are  quite  Involved,  whereas  numerical  calcula- 
max     ^  ' 

tlons  are  more  convenient  and  preferable. 

The  shock  thickness  Is  Infinitesimal,  and  It  is  always 
much  smaller  than  the  curvature  of  any  macroscopic  shock 
surface.   Consequently,  general  curved  MHD  shocks  In  steady 
three-dimensional  flows  can  always  be  approximated  locally 
by  a  planar  surface  for  which  the  present  solution  can  be 
used . 

The  solution  for  general  MHD  shocks  Includes  a  full  des- 
cription of  the  Interplay  of  several  characterizing  physical 
mechanisms;  Its  complicated  form  Is  hardly  surprising,  but 
solutions  for  many  different  limiting  cases  can  be  derived 
readily  from  it.   Final  quantitative  results  can  be  obtained 
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conveniently  only  with  the  aid  of  an  electronic  computer. 
In  this  regard,  a  subroutine  subprogram  (e.g.,  an  IBM  binary 
deck)  for  the  shock  solution  giving  values  of  all  downstream 
variables  with  the  values  of  upstream  parameters,  y,   and  the 
shock  strength  parameter  to  be  fed  in,  is  very  convenient 
in  solving  MHD  flow  problems,  while  a  table  containing  numeri- 
cal data  with  full  ranges  of  all  basic  parameters  covered  is 
not  feasible  here. 

After  a  complete  understanding  of  the  dynamical  stability 
criteria  for  general  MHD  shocks  is  acquired,  one  may  use  the 
present  shock  solution  with  assured  confidence  to  analyze 
quantitatively  many  MHD  flow  problems  which  occur  in  several 
branches  of  physics  of  great  current  interest,  e.g.,  the 
controlled  fusion  process,  geophysical  and  astrophysical 
problems . 


54.   A  table  containing  numerical  data  for  the  solution  of 

general  MHD  shocks  may  be  regarded  as  a  generalization  of 
the  familiar  table  for  the  solution  of  gasdynamic  shocks 
(cf.  NACA  Ames  Research  Staff,  Report  1135  (1953))-   How- 
ever, in  MHD  flows  the  presence  of  two  additional  indepen- 
dent parameters  due  to  the  magnetic  field  vector,  in 
addition  to  the  Mach  number,  renders  a  complete  table 
covering  full  ranges  of  all  characterizing  parameters  to 
be  entirely  unfeasible. 
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Appendix  A 

Using  (5a)  and  (5b)  to  eliminate  u  and  B  in  (2b,  c,  e, 
and  f)  and  expressing  them  in  terms  of  new  variables  we  get 

p_l     ^2  f  k  [(t-^)2(1+ts)2-(t-s)2(i+tj25)2]  ) 

^=-^  1-^+^ T^ 2 2 (Al) 

YMj    1+t   I     p  2M^T  (1+s  )(t-;2^)        J 

k  (/5-s)(1+t2)(t-s)     ^/l,^-p^ 

p[l  +  ^p p ]  =  llitlLl  (A2) 

MjT(1+s  ){^-j6)  T-e 

-  ^  T(;)-e)(T-s)  ^^3) 

s(t-£)(t-^) 

M^    .    k^s(l+Ts)    M?T^(l+e^)    ^    k-,(T-s)^(^-e)(l+T;25) 

p[^  +  ^  +  -^ 2—]  -  -^ ?—  +  -^  +  -^ 2 ?— 

2    Y-1     t(1+s  )      2p(T-£)       Y-1     (T-j25)  (T-e)(l+s  ) 

(A4) 

These  are  four  algebraic  equations  governing  t,  fi,    e,   V,    and  p, 
with  known  upstream  parameters  M_,  k  ,  s  and  the  given  gas 
adlabatlc  exponent  y.      e(T,p)  and  /!i(t,p)  may  be  readily  solved 
from  (A2)  and  (A3);  they  are  given  In  (8b)  and  (8c).   Elimina- 
tion of  j6   In  (Al)  by  the  use  of  (8c)  yields  P(t,p)  In  (8d). 
Using  (8b,  c,    and  d)  to  eliminate  e,   j6 ,    and  P  In  (a4),  we  get 

f(T,p)  =  0  In  (8a).   Elimination  of  e  In  (5a)  by  the  use  of 

-2    — 
(8b)  yields  u  (t^p)  In  (8e).   Using  (8c)  to  eliminate  j6   In  (5ti) 

yields  B^(t,p)  In  (8f ) . 
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Appendix  B 

Consider  a  quadratic  and  a  cubic  equation  of  x  as  follows: 

ax^  +  bx  +  c  =  0  (Bl) 

dx^  +  ex^  +  fx  +  g  =  0  (B2) 

Elimination  of  x  from  (Bl)  and  (B2)  yields  the  following 
equation  relating  the  coefficients: 

a[ (ag-ce)^  +  (af-be) (cf-bg) 

+  d[c-^d +bc(ag-ce)  +  (cf-bg)  (b^-2ac)]  =o       (B3) 

Since  (8a)  and  (8b)  are  cubic  and  quadratic  equations  of 
p  respectively,  p  can  be  easily  eliminated  by  the  use  of  (B3). 
The  derivations  are  simpler  If  we  use.  Instead  of  (8a)  and  (8b), 
the  equivalent  but  more  concise  expressions  In  (lla)  and  (lib). 
Taking  x  =  p/u  ,  we  determine  readily  the  expressions  of  the 
corresponding  coefficients  a,  . . . ,  g  by  comparing  (Bl)  and  (B2) 
with  (lib)  and  (lla)  respectively.   Substitution  of  these  Into 
(B3)  yields 

2/2   2  V  * 
b  bf(u  -b'^)f   =  0 
n  t  ^  n  n-* 

where 

*    *    u   u^   b   b, 
'     =  ^    (-'   -T-'   T-^   -T'   ^-'    y)    -    (^"n+-t)^^n[<-^(^-l)-^?] 

+  (eu^+u^)2(u^-eu^)jb^[ab^+(Y-i)u^b^][(Y-i)b2+(2-Y)u2]  + 
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+    (T+l)v2bJ  +  (u2-b2)|a[bJa+Yb2-(Y+l)b2]    +   y^^^,^b^ 


+  u^bjY\b2+  (^+l)(2-Y)bj2u^b^-u^bj]]| 

+    (Y+l)(u^-eu^)^iu^(u^u^-b^b^)2(2u^b^-u^bJ 

2      2  ? 

-    (u^-b^)[rY+l)u   b    (u   b,-u,b    )+(u  u,-b   b,  )(ab  +yu  "u.b. 


2  2  2     2 

with  a  E  2a  +  (y-1)u  .   Since  b  /^  0,  b,  ;^  0,  and  u  ;^  b  , 


in  general 


f*  =  0  .  (b4) 


It  gives  the  solution  for  compressive  MHD  shocks,  and  it  is 
a  polynomial  of  degree  nine  in  u  .   Expressing  (b4)  in  terms 
of  M_,  k  ,  s,  and  r  by  the  use  of  (lOa),  we  get  (17).   (b4) 
is  a  convenient  representation  of  the  shock  solution  for  one- 
dimensional  unsteady  flows. 


Appendix  C.   Normal  incidence  (u,  =  0) 

_o 
At  u,  =  0,  e,  fi,    and  u   in  (lib,  c,    and  e)  simplify  con- 

-   -      -2 
siderably  as  shown  in  Section  II,  whereas  p,  P,  and  B  in 

(lla,  d,  and  f)  remain  unchanged.   (B^)  also  reduces,  after 
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setting  u,    =   0  and  u     =  u   in  It ,    to 


u^I   =   0  (CI) 


where 


n  n         V  r       / 


+   e\^[2aV+  (y-l)uV][(i|-Y)u2+  (y-l  )b^]+(Y+l  )b  V 
-   Yu2b2J+   eb^f[2a2+(Y-l)u2][[2a2+(Y-l)u2-(Y+l)b2](u2-b2) 
+   b2[(2+Y)u2-Yb2]j  +   u2b2[27b2+(Y+l)[(2Y-l)b2+2(2-7)ij2]jj 
+    2(Y+l)b2b2[,;2^2_^^2_^2)^^2_^2)^     _ 

2 
Thus  1=0  Is  a  cubic  equation  of  u  . 


Appendix  D 

(a)   Aligned  flow 

B  =  const.  X  pu   in  aligned  flows,  and  hence  ii     =  q,  , 
ii^  ^  Q  ,    and  B  =  pu.   Ahead  of  the  shock  we  have 

2    2 
b    bf    ,2    k 
_n__t__b__l_  ,      . 

2     2  ""   2  ~  „  ?  ~  ^  '^^-^'' 

u^   "t   ^    ^I 

where  the  constant,  k,  is  defined.   (Dl)  can  be  used  to  elim- 
inate b^  and  b^  In  all  equations.   The  original  four  Independent 
parameters  (u^a,  u^/a,  b^a,  b^/a)  corresponding  to  (M,,  k  , 
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s,  t)  then  reduce,  in  this  special  case,  to  three  Independent 
parameters  (u^/a,  u^/a,  k)  corresponding  to  (m  ,  k  ,  t).   (lla) 


becomes 


Kp(p-l)u^  +  (l-Kp)i2pa^(l-Kp)-[Y+l-(Y-l)p](u^-Kpu^)   =  0  .    (D2) 
(llb)-(llf)  simplify  to 


e.6-.        """t'^"-'' 


2-2    -   2       ' 
^n     P^t  ~  ^P^ 


Y^^^(p-l)  J     Kpu^[l+  (1-2k)p]  ) 


pa^     I     2u^(l-Kp")' 
2 


Also  (b4)  simplifies  to 


(D3) 


f*(e,  ^;  -^,  k;  y)  =  K^u^[2a2+  {y-l)u^]A''   =  0  (d4) 

where 

A*  =  KTJ^(£U^+U^)^  KU^(£U^+U^)  +  U^  [  2 -7+ (  y-l  )  K  ]  (  U^-EU^  ) 

+  (1-k)(u  -£u,  )2U2a^+(Y-l)u^+K['^2_(^^3_)  2^1^^^  ^   . 


■X- 

A  =^0  gives  the  shock  solution,  and  it  corresponds  to  A  =  o 
in  (18). 
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(b)  Very  weak  magnetic  field 

Here  b  /a  «  1  and  b,/a  «  1,  whereas  b  /b.  Is  finite 
n  t  n   u 

We  get  from  (lla)-(llf) 

(t+1)u2       ^  ^^^u^{p-l)  _   P%-^t(V^t) 


(D5) 


2a^+(Y-l)u^  '     u^+pu^     •  ■    pu^+u^(byb^) 

T^n(p-^)    -2   ^n+P^^t    ^2    .  ,   P^-1 

P  =  1+^ ,   u  =   _2  2    .  B  =  1  +   ^   2/^2  • 

pa  p  u  1  +  byb^ 

-      -       -2 
Thus  In  this  limit  p,  e,   P,  and  u  are  unaffected  by  the 

presence  of  the  very  weak  magnetic  field;  they  are  Just  the 

well-known  result  for  a  gasdynamlc  shock.   Conversely,  the 

magnetic  field  behind  the  shock  being  characterized  by  0',  and 

-2 

B  ,  depends  on  the  fluid  properties  and  the  shock  strength 

parameter  p .   Since  P  >  1  for  real  shocks  satisfying  the 

-2 
entropy  condition,  we  have  B  >  1  always  In  this  limit. 

(b4)  becomes 
^    u   u,   b        b, 

^  ^^'  -i^'  -^'  -f«  1'  ^«  1'  '^)  (^^) 

^  "nK-^^t)'{^nt2a^+('V-l)u^]  +  u  Ju^b^-(2-y)u^b^  ]j  w"  =  0 
where 

^  O  Q 

W  =  [2a  +(Y-l)u^(eu^+u^)-(Y+l)u^u^(u^-eu^)  . 

e  does  not  enter  in  the  curly  bracket  of  (d6);  hence  the  shock 
solution  is  given  by  W  =0  which  corresponds  to  W  =  0  in  (19). 
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(c)   Very  strong  magnetic  field 

Here  b  /a  »  1  and  b  /a  »  1,  whereas  b  /b,  is  finite, 
(lla)  becomes 

pV|2paV  -  uV[y+1-(y-1)p]]  =   0  .  (D?) 

Apparently  p  =  o  is  not  a  physically  realistic  solution.   Then 
we  get  from  (D7)  for  p  /^  0  ,  together  with  (llb)-(llf), 

p-  -  (^+^)^n  ^  _   %(P-l)K-%(V"n)) 

2ay(l  +  b^/b^)  .  (Y-l)u2  '  '  "  u2+p-u2-u^u^(p--i)(b^/bj  ' 

,   ,   ,  2  -  2     (^^) 

^  ^  %(\/^n)   ;  -t  ^  ,  ^  ,  ^  TU^Zil  (,  ^  \) 


"n  +  ^^V^n)  pa  ^n' 

-2      1   f  2    -     -       ^t   /       -2 

"   ==  ^^-2r'n+^P"t-(p-^)%  — ]     '    B   =1. 
p  u  V  n   ^ 

— ?  — i    -^ 

Thus  we  have  /  =  s  and  B  =  1  which  means  B^  =  B^.   The  very 

strong  magnetic  field  remains  unchanged  across  the  shock  while 

-     -      -2 
the  fluid  variables  p,  e,    ?,    and  u   are  affected  by  the  presence 

of  the  magnetic  field  only  through  its  orientation  given  by 

b,/b  . 
V     n 

(b4)  becomes 

„    u   u,   b        b, 
f'    f  n    t    n^^T    t^^T    \ 

^  ^^'   —'   ^'   —  »  ^'   —»   1'  y)  (D9) 

=  b  [b  (eu  +u , )  +7b, (eu  +u, ) (u  -eu, )+(7+l)b  (u  -eu, )  IS  =  0 
n^  n^   n   t'   i^  t  ^   n   t'^  n    t '^  ^  '   ■'  n  ^  n    t'    ■' 

where 
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b^[2a^b^+  ('y-l)u^b2]  +  (y+l  )u^b2  (u^b^-u^b^  )  J  (u^-eu^  ) 

The  term  Inside  the  square  bracket  In  (D9)  corresponds  to 
p  =  Oj  so  It  Is  not  a  physical  solution.   Shock  solutions 
are  given  by  S   =0  which  corresponds  to  S  =  0  in  (20) . 
It  is  of  interest  to  observe  from  (d8)  that  p  and  P 

for  this  very  strong  magnetic  field  case  can  be  obtained 

2 
simply  by  replacing  a  in  their  corresponding  gasdynamic 

solution  with  a  /(l  +  b,/b  ). 


(d)   Zero  flow  deflection  angle 

e  =  0  in  this  case,  and  (lib)  gives 

u  (p-l)[u,  (u^-pb^)  +  pu  b  b,  ]  =  0   .  (DIO) 

n^^^'-t^n^n'^    '^nnt-'  ^    ' 

Let  us  first  consider  the  case  of  p  =  1 .   (Ha)  becomes 

2(u^-b^)z^  =  0  (Dll) 

where 

-7*    2,  2    ,2      2w  2  ,  2v 
Z-,  ==  u  b,  -  (u  -a  )  (u  -b  )  . 
1    n  t    ''  n    '^  ^  n  n'' 

2     2 
However,  e  is  not  necessarily  zero  if  both  p  =  1  and  u   =  b  , 

J  r  n    n' 

because  the  denominator  of  (lib)  becomes  zero  also  (see  IIB-l-g 
and-  Appendix  D-g  for  further  dlBcusslons ) .   We  consider 

here  only  u^  f   b^.   Then  Z-^  =  0  and  (llc)-(llf)  give  /5  =  s, 

-2    -2 
P  =  u  =  B  =1.   These,  together  with  e  =  0,  show  that  the 
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shocks  are  of  zero  strength  across  which  the  state  of  gas 

does  not  change  at  all.   Z-,  =  0  is  the  familiar  equation  for 

the  characteristics,  or  the  small  disturbance  wave  speeds 

(i.e.,  the  limiting  solution  of  shocks  of  very  weak  strength). 

Two  roots  of  u   in  Z-,  =  0  describe  fast  waves  whereas  the 

other  two  roots  describe  slow  waves  in  general.   These  are 

the  only  four  solutions  describing  shocks  of  zero  strength. 

Next  we  consider  the  case  that  the  square  bracket  in  (DIO) 

vanishes.   This  gives 

2 
u  u, 

"  ^  (D12) 


b  (u,b  -u  b,  ) 
n^  t  n   n  t  ^ 

which  describes  shocks  of  finite  strength.   Since  P  >  1  for  a 
physical  shock,  (D12)  yields  the  necessary  condition  that 

u  (u  u,  =  b  b,  )  >  u,b  (D13) 

n^nt    nt^—  tn  \      ^  i 

for  a  physical  shock.   The  corresponding  expressions  of  />,  P, 

-2      -2 

u  ,  and  B   can  be  easily  obtained  by  substituting  p  given  by 

(D12)  into  (llc)-(llf). 

Setting  e  =  0  in  (b4)  yields 

„    u   u,   b   b,  „  ^ 

f  (°-  ^^  ^'  ^'   ^■'    ^-J  =  2Z*Z*  =  0  (W4) 


where 


■X-  o      '7>      o  2  2 

Z-,  =  (u  b,-u,b  )[u,  (u  -b  )  -  yu  u.b  b,  -  (Y+l)u^b^  ] 
1    ^nt   tn'^t^n  n^   'ntnt^'   ■'nn-' 

2      2 
-  u  u,  [2a  b  +  (y-1)u  (u  b  +  u, b,  )  ]  . 
nf-     n   ^'   ^n^nn    tt^-' 
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Z*  =  0  corresponds  to  Z-j^  =  0  in  (21)  .   Zg  =  0  corresponds 
to  Z-|  =   in  (21).   In  the  limit  of  very  weak  magnetic  field 
(b  /a  «  1,  b^/a  «  1,  and  b^b^  finite),  Z^  =  0  gives  "^  =  a 
and  It  corresponds  to  sound  waves  In  gasdynamlcs,  while  Z-,  = 
gives  u,  =  0  (or  t  — >■  co) ;  hence  e  =  0  which  Is  the  solution 
for  a  normal  shock  wave  of  finite  strength. 


(e)   Normal  upstream  magnetic  field  (B-  .  =  0 ) 

Consider  this  first  as  the  limit  of  B,  ,  ^^  0  so  that 
e,  and  e   are  well  defined.   Setting  b,  — >•  0  In  (lla)  gives 

(u^-pb^)^  2pa2-u^[T+l-(T-l)p]   =  0  (D15) 

Let  us  consider  u  t^  pb  first.  The  vanishing  of  the  curly 
bracket  in  (D15)  yields  the  same  solution  of  p  as  that  for  a 

gasdynamlc  shock  given  in  (D5)-   Setting  b,  =  0  in  (llb)-(llf), 

—  1       —2  —       —2 

we  get  /5  =  -u,/u  =  -t   =  s,  B  =1,  whereas  z,    P,  and  u 
L*   n 

are  precisely  the  gasdynamlc  solutions  given  in  (D5)«  Thus 
we  have  for  this  case  a  pure  gasdynamlc  shock  with  B„  =  B^ . 
Obviously,  u,  remains  unchanged  across  this  shock  and  a  normal 

shock  exists  if  u,  =  0- 

2    -  2 
Next  we  consider  u  =  pb  which  is  a  double  root  here. 

n   "^  n 

Straightforward  substitution  of  b,  =  0  into  (llb)-(llf)  gives 
rise  to  a  set  of  indeterminant  equations.  A  precise  examina- 
tion of  the  limiting  process  is  necessary.   It  follows  from 

(lla)  that  as  b,  -^0  and  u^-pb^  -^>0 ,  b?(u^-pb^)"^  must 
t  n'^nt^n^n' 


-  68  - 


remain  finite.   To  approach  the  limit  properly  it  is  convenient 
to  introduce 

2-2 
b^  =  &x  and  ^^^-pb^^  =  &y  (Dl6) 

where  x  and  y  remain  finite  as  the  small  parameter  5  tends  to 

—  —       '7^    —  ^  —  — 

zero.   Denoting   ^  =  (b^b^p (p -1 ) )/(u^-pb^)  -  (xb^p (p-1 ) )/y, 

_      2/2 
we  get  from  (lla)  as  6  — >  0  (hence  b   -^ b  and  p  -^  u  /b  ), 

.2  ^\ 


^  =  -[(p"-i)[2(i-  ^)-{l-^){^-^)^]  (DiT) 

which  is  finite.   The  sign  in  front  of  the  radical  in  (Dl?) 
should  be  determined  by  the  physical  boundary  condition  in  the 
problem.   Using  (Dl6)  to  eliminate  b,  and  u   in  (llb)-(llf) 
and  then  taking  the  limit  of  6  — >0,  we  get 

u  [u,  (p-1)  +u^lj 


e  = 


"n  +  P"t+%"t^ 


1  +  ^   (p--l-  %^)    ,      u^   = 


a 


%^  -  "t 

%  +  "t^ 

( 

"n  +  (P"t  +  %^)^ 
-2    2 
P    u 

b2 

^-l^e, 

u _  +  u^^ 

(D18) 


where  u  =  pb  .   Since  u  /u,  =  t  =  -s   here,  ^  can  be  solved 
n   '^  n  t 

from  the  equation  for  ^   in  (Dl8)  as 


^  -  ^^^^-  =  ta.n{ii   -ii    )      or  ii  =   -\^  +tan~-^^  .  (D19) 

l+/5s        ^   ^ 

Hence  ^  represents  the  tangent  of  the  deflection  angle  of  the 

magnetic  field  across  the  shock.  Its  sign  can  be  easily 
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determined  here  as  required  by  the  boundary  conditions. 

2    2      2 
Since  B^.  =  B_   =  B_   ,  we  have  also  from  the  equation  for 

1      -i-  }Tl  r  ,  n 

B^  In  (Dl8) 

which  Is  equivalent  to  (D19)'   Thus  a  finite  tangential  mag- 
netic field  behind  the  shock  appears,  whereas  It  Is  absent 

ahead  of  the  shock;  this  type  of  shock  Is  called  a  complete 

^  6  13 
switch-on  shock  '    '        and  obviously  It  belongs  to  fast  shocks. 

Since  p  >  1  Is  required  by  the  entropy  Increase  condition 

2    -  2 
for  finite  shocks,  we  have  from  u  =  pb  that 

'  n   ^ 

u^  >  b^  .  (D21) 

Also  since  ^  >  0  is  necessary,  whereas  p  >  1  and  y  >  1 
always,  (Dl?)  yields 

b^  >  a^  +  ^  b^(p-l)  >  a^  .  (D22) 

(D21)  and  (D22)  are  two  necessary  conditions  for  the  existence 
of  a  complete  switch-on  shock.   Evidently  ^  Is  the  most  natural 
parameter  to  describe  a  switch-on  shock.   If  I  =0,  (D18) 
reduces  to  the  solution  describing  a  pure  gasdynamlc  shock 

as  discussed  previously. 

—2     —2       2       ? 

u^  =  (u  (I+te)  )/(l+e  )  by  definition  and  we  get,  after 

eliminating  e  by  the  use  of  e  in  (Dl8), 

|u   2 
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Hence,  normal  shock  does  not  exist  unless  ^  =  0  and  u,  =  o 

2    -  2 
In  terms  of  M.-^ ,    k^,  and  s,  u  =  pb   becomes 
II  n   "^ 

2 
M 
;r       I 


(D24) 


k^(l+S^) 

The  conditions  in  (D21)  and  (D22)  become 

M^  >  k^(l+s^)   and   k  >  1  ,  (D25) 

respectively.   Using  (D24)  to  eliminate  p  in  (D17)  yields 

r-[ ^  -  1]  2(1   1)-(y-1)[ ^-  -1]    (D26) 

k^(l+s2)      I     4        k^(l+s^)    J 

Substitution  of  (D26)  into  (D20)  yields  B^  .  =  Bp  ^(M, ,k  , s;y) 
in  (2^1). 

Setting  b,  =  0  and  b   ==  b  in  (b4)  yields 

^""^^^  -f'    "T'  I'  0'  '^)  "  ^V2  =  °  ^^^^) 

where 

N*  =  N*(e,  -^)  =[2a^+(Y-l)u^](£u^+u^)-(Y+l)u^u^(u^-£u^) 

■Jf      ^       v^       o  o 

N2  =  N2(e,  — )  =[b  (eu^+u^)  -  u^u^ (u^-eu^ ) ] 

+  (u^-eu^)2(u2-b2)[2a2+(Y-l)u^-(Y+l)b2]. 


•X- 


N-,  =  0  coincides  with  W  =  0  in  (d6)  for  a  gasdynamic 

2    —  2 
shock.   It  describes  the  case  of  u  7^  pb   and  corresponds  to 

*  2—2 

N-|  =  0  in  (22).   N^  ^   0  describes  the  case  of  u  =  pb   for 
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switch-on  shocks.   It  can  also  be  obtained  readily  by  elim- 
inating ^  In  the  equations  for  e  In  (Dl8)  with  the  aid  of 
(DI7)  and  setting  p =  u  /b  .  N^   =  0    corresponds  to  Ng  =  0 
In  (22). 

For  the  special  case  of  aligned  flows  that  u,  =0,  e, 
^,   and  u   simplify  to 

e  =  j6  ^  ^        and   u^  =  ^^  ,  (D28) 

P 

-      -2 
whereas  P  and  B  remain  unchanged  as  given  In  (DI8) .   These 

agree  with  the  result  obtained  by  Frledrlchs,  and  Bazer  and 

Erlcson. 

If  B_  X  e   =0  strictly,  there  Is  no  unique  way  of 

defining  e,  and  also  e  .   A  degenerate  two-parameter  family 

of  solutions  exists.   For  the  case  of  a  pure  gasdynamlc  shock 

(N-,  =  0  j ,  the  determination  of  e,  and  e   are  entirely  Irrelevant 

For  the  case  of  a  complete  switch-on  shock  (Np  =0),  the 

finite  magnetic  field  component  on  the  transverse  plane 

,  -X   ,  _A 

(normal  to  e  }  downstream,  B„  „?  "lay  be  directed  at  any  angle 

n  •t'  ^  i 

on  the  boundary  conditions.   Then  we  have  e,  I |b„  „  (=  B„  , ), 

and  e   Is  also  defined. 
z 


A  switch-on  shock  differs  from  an  Alfven  shock  In  that 

2 
'T 


p  7^  1  and  B„  j^  1 .   It  Is  a  compressive  shock  and  entropy 


2    -  2 
increases  across  it.   However,  u  =  pb   means  simply 

2      2 
"p  n  ^  ^F  n'  '^^'^^   3-  switch-on  shock  may  be  followed  immed- 
iately by  a  non-compresslve  Alfven  shock  across  which  B„  ,  may 
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further  rotate  through  an  arbitrary  angle  in  the  transverse 
plane,  i.e.,  the  (t,z)  plane,  with  its  strength  unchanged 
while  the  flow  velocity  changes  both  its  magnitude  and  direc- 
tion accordingly. 


(f)   Tangential  magnetic  field  (B-^   =  q  ) 

_L  ,n 

Setting  b   =0,  hence  b,  =  b  in  (lla)-(llf)  yields 

2pa^  -  u^[y+1-(y-1)p]  +  pb^[Y+(2-Y)p]  =  0 

""n^^P"^)     ^     ,   X     -      T^n^P-l)  ,    b^P(P+l)   ,    , 

e  =      2-2         '      ^   "  ^^^^^    '      P  =  1  +  -^ fl 5 ^  ^^29) 

u^  +  pu^  a  p         2u^ 

2    -  2 

-2    ^n  +  P^t      =2        -2 
u   =  — 3^-2 —  ,  B  =  p   . 

p  u 


The  first  of  the  above  equations  is  a  quadratic  equation 

—  — ?     — ?      P      ? 

of  p  now.   Since  u,  =  (u  (I+te)  )/(l+e  )  by  definition,  we 

get,  by  the  use  of  (D29),  u,  =  1  always,  as  expected.   Normal 

shocks  exist  at  u,  =  0  for  which  (-D29,)  yields  e  =  o  and 

-2    -2    1 

u  =  u  =  -^  with  the  other  expressions  unchanged. 

P 
Setting  b   =  0  in  (B^)  yields 

^*^^'   X-'  -t'    0'  ■^'  ^^    ^    (2-7)u^b(u^-eu^)T*  ^   q       (D^O) 
where 

T*  =  (Y-2)u^b^(£U^+u^)^-u^[2a2+(Y-l)u^+Yb^](u^-eu^)(eu^+u^) 
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it 

T     =0    corresponds  to  T  =  0    in    (26). 


(g)   Non-characteristic  solutions  with  p  =  1 

It  was  shown  In  Appendix  D-d  that  four  solutions  of  t 
at  p  =  1  correspond  to  shocks  of  zero  strength  across  which 
the  state  of  gas  does  not  change  at  all.   They  give  the 
characteristic  solutions  and  may  be  regarded  as  the  limiting 
solutions  of  very  weak  shocks.   However,  f (t,p)  =  0  as  given 
in  (8a)  is  a  polynomial  of  degree  six  in  t.   For  any  given  p, 
six  real  and  distinct  roots  of  t  exist  in  general.   The  other 

two  solutions  of  t  at  p  =  1  are  given  by  taking  the  limit  of 

-  2      2 

p  —^1  and  u  — >b   simultaneously,  whereas  the  latter  condi- 

r      '^  n     n 

tion  has,  so  far,  been  assumed  to  be  not  valid. 

To  obtain  solutions  for  this  case,  it  is  convenient  to 

2    2 
define  p  =  l+5x,  u  =b  +6y  where  x  and  y  remain  finite 

as  6  — >0.   Substituting  these  into  (lla)  and  taking  the  limit 

of  6  — ^  0,  we  get 

y  =  I  b^x  .  (D31) 

Since  5x  >  0  is  necessary  according  to  the  entropy  condition, 
u^  -  pb^  -^  -5x(b^/2)  <  0  in  the  limit.   Thus  only  slow  shocks 
are  associated  with  this  case.   Similarly  we  get  from  (llb)- 
(llf),  with  the  aid  of  (D3l),  the  following  solutions  in  the 

limit: 

2 
u,      b       u,    -I 

^  =  [^-5-57(1  +  4)]  (W2) 

n      t      u 

n 

-  7^  - 


"u  b,  +  u,  b 

u,  b,  -  u  b 
t  t    n  n 

b   2 

u  n 

P  =  B^  =  1  .  (D35) 

Unless  b,  =  0  which  corresponds  to  a  special  case,  e  ^      , 

p  ^   s,    and  u  /=  1  in  general.   Hence  u„  7^  u-  and  B„  ;^  B.^..   The 

-    -    -2 
shock  Is  of  finite  strength  even  though  p=P=B  =1.   It 

is  not  appropriate  to  consider  p  as  a  shock  strength  parameter 

here  for  it  would  easily  lead  to  the  misleading  conception 

that  p  — >1  means  shocks  of  zero  strength. 

a  does  not  appear  in  (D31)-(D35)-   Therefore  the  Inertial 

energy  is  irrelevant.   The  change  of  state  across  this  type 

of  shock  is  essentially  kinematic  in  nature,  f)   in  (D33)  means 

ii^  =   2p-/r-^^  .  (D36) 

_2 
This,  together  with  B  =1,  shows  that  across  this  shock  the 

change  of  magnetic  field  occurs  only  by  reversing  its  tangential 

component,  i.e.,  B„  ,  =  -B-r  ,  .   It  has  precisely  the  same  effect 

p  ,  t     J-  J  ^ 

as  the  change  across  an  Alfven  (or  transverse)  shock  with  the 
tangential  magnetic  field  being  rotated  through  an  angle,  tt, 
in  the  transverse  plane,  together  with  the  accompanying  change 
of  the  flow  velocity  accordingly. 

Expressing  u  =  b   and  (D32)  in  terms  of  t,  M^,  k  ,  and  s, 
we  get  (28a)  and  (28b). 
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It  Is  also  of  interest  to  see  that  for  aligned  flows, 
(D32)-(D35)  reduce  to 

2u  u, 
^=  t^  =   -^Jl4  (D57) 

"t  -% 
P  =  u2  ^  B^  =  1  (D38) 

simply;  they  are  called  improper  slow  shocks  in  ref.  6.   The 
only  changes  across  the  shock  in  this  case  are  the  reversal 
of  the  tangential  components  of  both  flow  velocity  and  magnetic 
field. 
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Figure  Captions 

Pig.  1:   Velocity  and  magnetic  field  vectors  on  both  sides 
of  a  general  oblique  MHD  shock. 

Pig.  2:   The  shock  adlabatlcs  on  {^  ,Q)    plane  for  y  =  5/5, 
M^  =  3.0,  and  i^     =   ^0    .   The  values  on  the  upper 
right  corner  of  each  graph  signify  those  of  (M  ,  k  , 
1^  ;  y)  respectively,   co  and  co_  denote  the  fast  and 
slow  characteristics  (or  waves)  respectively.   r 
denotes  the  non-characterlstlc  solution  with  p  =  1. 

Pig.  3:   Pig.  (3a)  shows  the  positions  of  the  terminals  of  -u" 
relative  to  the  characteristic  locus  at  k  =  0.5 
and  2.0,  respectively,  for  each  of  the  cases  shown  In 
Pigs.  (3b-n).   All  speeds  are  normalized  with  respect 
to  a.     there.   Pigs.  (3b-n)  Illustrate  the  shock  adla- 
batlcs on  the  (P,0)  plane  for  y  ~   5/3,  k  =  0.5  and 
2.0,  ify  -   5  and  0  ;  the  values  on  the  upper  right 
corner  of  each  of  them  signify  those  of  (M-j-,  k-p, 
ij/^;    y)  respectively,   o)  and  co_  denote  the  fast  and 
slow  characteristics  (or  waves)  respectively.   r  denotes 
the  non-characterlstlc  solution  with  p  =  1. 
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This  report  was  prepared  as  an  account  of 
Qovernment  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  Implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  Information  contained  In  this  report, 
or  that  the  use  of  any  Information, 
apparatus,  method,  or  process  disclosed 
In  this  report  may  not  Infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  Information,  apparatus, 
method,  or  process  disclosed  In  this 
report. 

As  used  In  the  above,  "person  acting  on  behalf 
of  the  Commission"  Includes  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  Infor- 
mation pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with 
such  contractor. 
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